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ABSTRACT 

The  thesis  compares  the  analytical  solution,  two  marine  classification  society  design 
rules,  and  numerical  analysis  against  experimental  results  for  predicting  the  failure  modes 
(general  instability,  axisymmetric  buckling,  and  asymmetric  collapse  of  the  shell)  and  failure 
pressures  of  ring-stiffened  cylindrical  shells. 

The  analytical  solution  is  first  summarized  based  on  several  sources.  Design  rules  for  the 
classification  societies  are  then  presented  with  brief  explanations  for  each  one.  The  design  rules 
used  are:  American  Petroleum  Institute  ( Bulletin  on  Stability  Design  of  Cylindrical  Shells,  API 
Bulletin  2U,  Second  Edition,  October  2000)  and  Det  Norske  Veritas  ( Buckling  Strength  of  Shells, 
October  2002).  The  numerical  analysis  was  performed  using  the  software  package,  Method  For 
Analysis  Evaluation  and  Structural  Optimization  (MAESTRO™,  version  8.5,  Proteus 
Engineering). 

The  United  States  Navy  Naval  Sea  Systems  Command,  Submarine  Structural  Integrity 
Division  supplied  experimental  data  for  four  test  cylinders  that  covered  the  failure  modes  and 
allowed  comparison  between  experimental  and  analytical  /  numerical  results. 

The  comparison  of  experimental  to  predicted  data  found  the  design  rules  and  numerical 
solution  performed  adequately  in  predicting  asymmetric  buckling  and  general  instability  failure 
modes,  but  the  predictions  for  failure  pressure  were  unsatisfactory.  The  design  rules  were  overly 
conservative  in  their  predictions  of  failure  pressure  due  to  the  semi-empirical  solutions  used  in 
the  rules.  The  numerical  solution  was  only  slightly  better  for  the  same  failure  pressure 
predictions.  The  results  indicate  the  predicted  failure  pressure  for  a  cylinder  is  closely  tied  to  the 
size  and  dimensions  of  the  cylinders  used  for  determining  the  empirical  solutions.  These  results 
should  be  further  explored  to  determine  causes  and  corrections. 
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CHAPTER  1:  Introduction  and  Problem  Statement 

Recent  interest  in  submersibles,  submarines,  and  off-shore  drilling  rigs,  has  led  to  an 
increasing  demand  for  structural  design  information  on  ring-stiffened  cylindrical  shells  under 
uniform  external  pressure.  The  submarine  designer  today,  has  many  analytical  tools  and 
methods  available  to  help  determine  an  optimum  design.  The  widespread  use  of  ring-stiffened 
cylinders  in  the  marine  industry  has  resulted  in  a  significant  amount  of  interest  and  activity  being 
devoted  to  determining  the  failure  pressure  and  characteristics  of  these  cylinders.  Marine 
Classification  Societies,  such  as  the  American  Petroleum  Institute  (API),  the  American  Bureau  of 
Shipping  (ABS),  Det  Norske  Veritas  (DNV),  Germanischer  Lloyd  (GL)  and  others  have 
promulgated  design  rules  to  provide  guidelines  on  designing  and  building  stiffened  cylinders  for 
marine  use.  Other  research  has  been  conducted  using  numerical  methods,  such  as  finite  element 
analysis,  to  design  and  validate  the  structural  adequacy  of  these  ring-stiffened  cylinders.  By 
comparing  these  design  methods,  classification  society  design  rules  and  numerical  methods,  with 
experimental  results,  the  submarine  designer  can  obtain  a  better  understanding  of  the  strengths 
and  limitations  of  each  method. 

1.1  Definition  of  Failure  Modes 

Any  discussion  of  cylinder  failure  analysis  must  first  include  definitions  of  the  different 
failure  modes.  There  are  primarily  three  failure  modes  for  ring-stiffened  cylinders.  They  are 
axisymmetric  yielding  (AX)  of  the  shell  between  stiffeners,  asymmetric  buckling  of  the  shell 
between  stiffeners  (Lobar  buckling)  (L),  and  general  instability  (GI)  of  the  shell  and  stiffeners. 
Axisymmetric  yield  is  characterized  by  an  accordion  type  pleat  extending  around  the  periphery 
of  the  cylinder,  and  generally  occurs  when  the  shell  is  relatively  heavy  and  the  frames  are  closely 


spaced.  Lobar  buckling  is  characterized  by  inward  and  outward  lobes  or  dimples,  which  may  or 
may  not  develop  around  the  entire  periphery,  and  normally  occurs  when  the  shell  is  relatively 
thin  and  the  frames  are  strong  and  widely  spaced.  General  instability  is  characterized  by  the 
failure  of  both  the  shell  and  ring  frames  resulting  in  a  dished-in  surface.  General  instability 
normally  occurs  when  the  cylinder  is  relatively  long,  the  shell  is  thin,  and  the  frames  are  light. 

1.2  Literature  Search 

The  failure  of  cylinders  exposed  to  external  pressure  has  been  studied  for  over  a  hundred 
years.  As  early  as  the  1850’s,  attempts  were  made  to  understand  cylinder  behavior  by  using 
experiments  and  empirical  relationships  [1].  The  first  analytic  solution  for  a  non-reinforced 
cylinder  was  presented  by  G.  H.  Bryan  in  1888  [2].  During  this  time  period,  non-reinforced  flues 
were  observed  to  fail  in  fire-tube  boilers  at  a  pressure  much  less  than  the  hoop  stress,  which  led 
to  a  significant  amount  of  research  and  interest  in  the  subject.  As  a  solution  to  this  problem, 
stiffening  rings  or  bulkheads  were  added  to  reduce  the  unsupported  length  of  the  tube  [3].  The 
first  analysis  of  a  reinforced  cylinder  appeared  in  1913  by  R.  V.  Southwell,  followed  a  year  later 
by  a  solution  to  the  elastic  buckling  of  a  thin  shell  proposed  by  von  Mises  [1].  In  1934 
Widenburg  proposed  a  solution  for  asymmetric  buckling  that  was  independent  of  the  number  of 
lobes  of  failure,  which  made  the  solution  easier  to  calculate  [1].  Solutions  for  axisymmetric 
yield  were  first  put  forward  by  von  Sanden  and  Gunther  in  1920  [2].  In  1930,  Viterbo  presented 
a  modified  version  of  Sanden  and  Gunther’s  solution  [2].  Finally,  Pulos  and  Salerno 
incorporated  the  previous  work  and  presented  a  solution  that  included  the  Sanden  and  Gunther 
solution,  the  Viterbo  modification  and  a  term  to  account  for  the  bending  stress  in  the  cylinder 
caused  by  the  axial  pressure  [2].  For  elastic  general  instability,  the  first  reported  analysis  was 


presented  by  Tokugawa  in  1929.  In  1954  A.  R.  Bryant  developed  a  similar  equation  using  a 
different  methodology  [1]. 

Analytical  work  from  the  1950’s  onward  has  focused  on  obtaining  solutions  for  different 
boundary  conditions  and  more  fully  reconciling  the  analytic  predictions  with  experimental 
results  and  more  fully  understanding  the  effects  of  initial  cylinder  imperfections.  With  the 
advent  of  the  digital  computer,  programs  like  BOSOR  5  were  developed  that  could  use 
numerical  solutions  to  quickly  and  accurately  predict  failure  pressures  [1].  Further  developments 
relating  to  numerical  solutions  led  to  the  design  of  finite  element  programs,  like  AB  AQUAS™, 
that  could  provide  accurate  stress  and  strain  values  for  analyzing  cylinder  designs  [1]. 

1.3  Previous  Work 

Tighter  budgets  in  both  industry  and  government  have  forced  many  large  organizations  to 
look  for  cost  saving  measures.  One  such  perceived  cost  saving  measure  has  been  the  outsourcing 
of  many  functions  that  were  previously  done  within  an  organization.  An  example  of  this  is 
found  in  the  greater  role  that  marine  classification  societies  are  playing  in  certifying  and 
classifying  naval  vessels,  not  only  for  commercial  interests,  but  also  for  governments.  This 
interest  has  led  many  classification  societies  to  develop  extensive  rules  for  certifying  naval 
vessels  and  other  marine  structures.  These  rules  can  also  be  valuable  tools  for  the  submarine 
designer. 

In  a  recent  review  of  these  classification  society  rules,  D.J.  Price  used  two  marine 
classification  design  rules  and  compared  them  with  analytical  and  experimental  results  for  ring- 
stiffened  cylinders  [4].  His  work  indicated  that  the  two  rules  used  (ABS  and  GL)  were  accurate 
for  predicting  axisymmetric  yielding  and  lobar  buckling  when  compared  to  experimental  results. 


However,  they  did  not  accurately  predict  failure  by  general  instability.  Further  study  was 
indicated  in  this  area. 

1.4  Problem  Statement 

In  today’s  fiscally  constrained  environment,  the  submarine  designer  is  faced  with  the 
challenge  of  providing  the  best  structural  design  possible  at  the  lowest  cost.  Detailed 
confirmation  models  can  increase  costs  not  only  through  expensive  fabrication  but  also  time 
delays  for  constructing  and  testing  the  models.  If  the  designer  can  use  some  of  the  tools 
available,  like  classification  society  design  rules  and  numerical  solutions,  to  reduce  or  eliminate 
some  of  the  confirmation  models,  there  are  significant  cost  savings  to  be  anticipated. 

This  thesis  used  three  of  the  design  tools  available  (classification  society  design  rules, 
numerical  analysis  tools,  closed-form  analytic  solutions)  to  determine  the  failure  modes  and 
pressures  for  four  experimentally  tested  ring-stiffened  cylinders.  The  results  from  the  design 
tools  and  the  experiments  were  compared  to  determine  die  applicability  and  usefulness  of  these 
tools. 

This  thesis  was  not  an  exhaustive  study  of  classification  rales  or  of  numerical  analysis 
tools,  rather  it  was  an  application  of  the  design  tools  available.  Comparisons  and  conclusions 
were  drawn  based  on  the  results  in  order  to  provide  the  submarine  designer  a  better 
understanding  of  the  limitations  of  each  design  method. 


CHAPTER  2:  Approach 


For  this  thesis,  emphasis  was  placed  on  exploring  how  various  classification  society 
design  rules  predicted  failure  of  cylinders  that  were  similar  in  design  to  modem  submarine  hulls. 
Similar  design  meant  that  the  shell  was  relatively  thick  compared  to  the  diameter  of  the  cylinder. 
For  comparison  purposes,  a  numerical  analysis  was  also  performed  on  the  same  cylinders  using  a 
numerical  analysis  tool.  In  order  to  compare  results  with  previous  work,  experimental  failure 
data  was  obtained  from  the  Naval  Sea  Systems  Command  (NAVSEA)  Submarine  Structural 
Integrity  Division  on  the  same  test  cylinders  used  in  [4].  For  consistency  of  analysis,  the  scope 
was  limited  to  examining  ring-stiffened  cylindrical  shells.  The  test  cylinders  used  were  selected 
to  cover  all  three  modes  of  cylinder  failure,  allowing  for  comparison  of  not  only  failure  pressure 
but  also  failure  mode. 

2.1  Analysis  Techniques 

2.1.1  Analytical  Methods 

For  the  purpose  of  this  thesis,  the  analytical  methods  include  the  classification  society 
design  mles  and  the  closed-form  analytic  solutions.  These  analytical  methods  were  programmed 
into  MATHCAD™  for  consistency  of  approach,  clarity  of  symbolic  representation,  and  ease  of 
calculation.  Dimensions  were  input  into  each  computer  code,  which  provided  failure  pressures 
for  each  mode  of  failure.  The  lowest  calculated  pressure  was  considered  the  failure  pressure 
with  a  corresponding  failure  mode.  The  failure  modes  and  pressures  were  compared  to 
experimental  results  with  primary  emphasis  being  placed  on  agreement  of  failure  mode  and 
secondary  emphasis  on  failure  pressure. 


2.1.2  Numerical  Method 


To  determine  a  numerical  solution  for  the  failure  mode  and  pressure  of  a  ring-stiffened 
cylinder,  a  numerical  analysis  tool  was  used.  Analysis  was  performed  using  the  Method  for 
Analysis  Evaluation  and  Structural  Optimization  (MAESTRO™),  version  8.5  distributed  by 
Proteus  Engineering.  Models  of  the  test  cylinders  were  created  in  MAESTRO™  and  subjected 
to  increasing  submergence  pressure  until  failure  occurred.  The  associated  failure  mode  and 
pressure  were  considered  the  failure  point  for  the  model.  Once  again,  failure  modes  and 
pressures  were  compared  to  experimental  results  with  primary  emphasis  being  placed  on 
agreement  of  failure  mode  and  secondary  emphasis  on  failure  pressure. 

2.2  Design  Rules  Examined 

There  were  two  classification  society  design  rules  examined:  The  American  Petroleum 
Institute  (Bulletin  on  Stability  Design  of  Cylindrical  Shells,  API  (Bull  2U),  Second  Edition, 
October  2000)  [5]  and  Det  Norske  Veritas  (Recommended  Practice  on  Buckling  Strength  of 
Shells,  DNV-RP-C202,  October  2002)  [6].  The  specific  classification  societies  were  selected 
due  to  their  widespread  use  throughout  the  world  and  the  availability  of  documented  rules  for 
ring-stiffened  cylinders.  Additionally,  API  was  selected  because  of  its  widespread  us  in  the  U.S. 
while  DNV  was  selected  because  of  its  widespread  use  in  Europe.  By  using  these  two 
classification  societies,  a  concise  snapshot  of  guidance  relating  to  cylinder  design  could  be 
obtained  for  a  large  segment  of  the  marine  industry. 


CHAPTER  3:  Basics  of  Ring-Stiffened  Cylindrical  Shells 

The  main  structural  body  of  most  submarines  and  submersibles  today,  is  constructed  of  a 
cylindrical  parallel  mid-body  section.  These  cylinders  are  reinforced  with  ring-stiffeners 
(frames)  to  provide  additional  strength  to  the  shell  that  would  collapse  very  easily  if  not 
reinforced.  A  strong  cylindrical  structure  is  required  for  the  large  pressure  differential  between 
external  hydrostatic  pressure  and  internal  pressure  (normally  maintained  close  to  atmospheric 
pressure). 

3.1  Nomenclature 

Each  of  the  analytical  and  numerical  methods  incorporated  in  this  study  used  slightly 
different  terminology  for  cylinder  geometries  and  properties.  When  the  analytical  methods  were 
programmed  into  MATHCAD™,  the  symbols  used  by  the  source  document  were  generally  used 
in  the  program  to  avoid  confusion  between  the  published  classification  society  rules  and  the 
programs.  All  of  the  analytical  methods  required  the  calculation  of  the  moment  of  inertia  of  a 
combined  plate  and  stiffener  (Ie)  using  an  effective  shell  length  (Le).  The  formula  for  Ie  came 
from  [7],  while  the  formulas  for  Le  were  normally  contained  within  the  classification  society 
rules.  All  stresses  and  pressures  are  in  pounds  per  square  inch  (psi),  lengths  are  in  inches  (in), 
areas  are  in  square  inches  (in2)  and  moments  of  inertia  are  in  inches  to  the  fourth  (in4). 

3.2  Dimensions 

The  dimensions  of  interest  for  analyzing  ring-stiffened  cylinders  are  related  to  the 
cylinder  (shell)  itself  and  the  ring-stiffeners  (frames).  Terms  and  definitions  are  listed  below  and 
represented  in  Figure  1  and  Figure  2. 


1)  Cylinder  Length  (Lb):  Overall  length  of  the  cylinder  between  supports. 

2)  Radius  of  Cylinder  (R):  Mean  shell  radius. 

3)  Shell  Thickness  (t,  h):  Thickness  of  cylinder  shell  plating. 

4)  Ring  Spacing  (Lf  or  Lr):  Unsupported  length  of  shell  from  centerline  to  centerline  of  frames. 

5)  Web  Height  (hw):  Length  of  the  web  from  shell  to  the  shell  side  of  flange. 

6)  Web  Thickness  (tw):  Thickness  measured  across  web. 

7)  Flange  Breadth  (bf):  Width  of  the  flange. 

8)  Flange  Thickness  (tf):  Thickness  of  flange  measured  petpendicular  to  breadth. 

9)  Faying  Width  (b):  Contact  width  of  ring  frame  to  shell,  normally  equal  to  tw. 

10)  Effective  Shell  Length  (Le):  Usually  some  fraction  of  Lf  specified  in  the  individual  solution. 

1 1)  Area  of  Stiffener  (Af  or  Ar):  Cross-sectional  area  of  the  ring-stiffener. 

12)  Effective  Area  of  Stiffener  (Aetr  or  AT):  Cross-sectional  area  of  combined  stiffener  and  Le 
of  shell. 

Figure  1:  General  Cylinder  Dimensions 
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Figure  2:  General  Stiffener  Dimensions 


3.3  Stresses  in  Cylinders 

Stresses  in  cylindrical  pressure  vessels  must  be  discussed  briefly  in  order  to  provide  a 
background  for  the  derivation  of  the  analytical  solution.  To  begin  with,  a  cylinder  can  be 
considered  a  thin-walled  shell  if  the  ratio  of  the  radius,  R,  to  shell  thickness,  t,  is  greater  than  ten. 
With  this  assumption,  the  determination  of  the  stresses  can  be  accomplished  using  statics  alone. 
All  of  the  cylinders  under  consideration  for  this  thesis  are  treated  as  shells.  Another  assumption 
in  the  analysis  is  that  hydrostatic  pressure  is  considered  constant  across  die  shell. 

From  classic  static  analysis  it  can  be  shown  that  cylindrical  shells,  exposed  to  hydrostatic 
pressure,  have  two  basic  stresses  imparted  to  them  by  the  pressure:  hoop  stress  and  axial  stress 
[8].  The  equations  for  these  stresses  are  shown  below: 


1)  Hoop  Stress: 

3 

ii 

0) 

2)  Axial  Stress: 

21 

(2) 

Where  p  is  defined  as  the  external  (or  internal)  pressure,  R  is  the  mean  shell  radius  and  t  is  the 
thickness  of  the  shell. 

Once  the  shell  is  stiffened  using  ring-frames,  the  hoop  stress  analysis  becomes 
complicated  because  non-uniform  deformation  of  the  shell  is  introduced  in  the  radial  direction. 
Additionally,  there  is  a  beam-column  effect  due  to  the  pressure  acting  in  the  axial  direction.  The 
effects  introduced  by  adding  ring-frames  are  discussed  in  detail  in  Chapter  4. 


<• 
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CHAPTER  4:  Analytic  Solution 

While  all  failure  modes  are  addressed  individually,  there  was  no  comprehensive 
theoretical  solution  that  addressed  all  modes.  Reference  [1]  provides  a  good  summary  of  the 
current  closed-form  analytic  solutions  that  are  widely  used. 

When  trying  to  determine  how  a  cylinder  will  fail,  it  is  often  advantageous  to  look  at 
some  key  parameters.  A  first  indicator  of  the  failure  mode  of  a  cylinder  is  found  by  plotting  the 
cylinder’s  slenderness  ratio  (A)  against  the  pressure  factor  (vp)  [9].  A  has  the  following 


nondimensional  value. 


3  u; 


\|/  is  the  ratio  of  the  shell  buckling  pressure  (pc)  to  the  hoop  pressure  at  yield  (py). 


Pc 

V  =  — 


For  most  steel  cylinders,  the  following  assumptions  can  be  made;  v=0.3  and  L/2R  » t/2R.  By 
making  these  assumptions  the  equation  for  vp  becomes  [9]: 


A  plot  of  \| /  verses  A  is  shown  in  Figure  3. 
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Figure  3:  Failure  Pressure  Ratio  versus  Slenderness  Ratio 


If  the  slenderness  ratio  is  less  than  roughly  1.14  then  the  cylinder  should  fail  by  axisymmetric 
yield  (AX),  and  when  it  is  greater  than  1.14  it  should  fail  by  lobar  buckling  (L).  If  the  shell  and 
stiffeners  are  not  sufficiently  sized,  the  cylinder  may  fail  by  general  instability  at  a  pressure  less 
than  that  predicted  by  the  \\ /  verses  A  curve. 

Another  very  important  factor  for  the  analytic  solutions  is  the  treatment  of  boundary 
conditions.  The  literature  devotes  a  significant  amount  of  research  and  discussion  on  what  types 
of  boundary  conditions  to  use  for  analysis,  with  methods  ranging  from  fully  clamped  to  simply 
supported  ends.  In  reality,  both  extremes  are  difficult  to  create,  so  the  experimental  results  fall  in 
a  range  between  the  two  extremes.  For  this  thesis,  no  discrete  boundary  conditions  were 
required  as  inputs  to  the  equations  because  the  analytic  solutions  used  do  not  distinguish  between 
differing  boundary  conditions. 

4.1  Axisymmetric  Yield 

Axisymmetric  yield  has  been  studied  since  the  1920’s.  As  discussed  in  Chapter  1.2, 

Pulos  and  Salerno  presented  a  closed-form  solution  for  axisymmetric  yield  in  1961.  It 


♦ 
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incorporated  previous  works  of  van  Sunden  and  Gunther,  and  Viterbo  and  includes  a  previously 
neglected  beam-column  effect  due  to  hydrostatic  pressure  acting  in  the  axial  direction  of  the 
cylinder  [10].  The  governing  differential  equation  for  the  Pulos  and  Slemo  equation  is: 


j4  i 


(T  '  p-R  d2  E-t  _  f.  v] 

7W  +  - - :W+ - W  -  p-  1 - 

v«4  ;  >  *!  E2  l 


(6) 


Where  w  is  the  radial  displacement  and  D  is  the  flexural  rigidity  of  the  shell  and  is  defined: 


D:= 


Et 


12-(l  -  v2) 


(7) 


The  beam-column  effect  term  is  —  which  makes  equation  (6)  a  non-linear  function  of 

2 

pressure.  This  term  was  neglected  in  the  previous  analyses  of  axisymmetric  yield  and  greatly 
improved  the  accuracy  of  the  results.  For  deriving  the  governing  equations,  a  coordinate  system 
for  a  shell  element  is  used  in  reference  [10]  and  is  shown  in  Figure  4. 


Figure  4:  Element  of  a  Cylindrical  Shell 


In  order  to  solve  the  non-homogeneous  differential  equation,  the  general  solution  of  the 
governing  equation  was  written  as  the  sum  of  the  solution  of  the  homogeneous  equation  and  a 
particular  solution  [1].  The  solution  to  the  homogeneous  equation  produces  four  roots  (Xi,  X2, 
X3,  X4).  By  analysis,  placement  of  the  origin  of  the  coordinate  system  to  take  advantage  of 
symmetry,  and  trigonometric  identities,  the  general  solution  can  be  written  as: 

w  =  B  cosh  +  F  cosh  (l-~)  (8) 


where  B  and  F  are  new  arbitrary  constants  of  integration  [10].  After  further  mathematical 
substitutions,  several  dimensionless  parameters  were  introduced  into  the  solution  to  allow  ease  of 
solving  the  problem.  Four  of  these  dimensionless  parameters  (Fi,  F2,  F3,  F4)  were  transcendental 
functions  based  on  the  geometry  of  the  cylinder.  Pulos  and  Salerno  graphed  these  transcendental 
functions  in  reference  [10]  to  allow  a  quick  solution  to  be  found  for  a  cylinder  with  known 
dimensions.  Finally,  an  equation  for  the  failure  pressure  of  the  cylinder  was  determined  The 
Poulos  and  Salerno  equation  is  used  in  this  analysis  and  is  shown  below: 


Where: 
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a  := 


a  +  P  +  (l  —  p)-F  ! 
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Once  the  variables  were  defined,  an  iterative  process  was  required  for  the  general  case  where  the 
parameter  y  was  not  zero.  Iteration  was  begun  by  assuming  y  was  zero,  and  then  finding  the 
corresponding  failure  pressure.  Having  this  interim  failure  pressure,  y  was  recalculated  solving 
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the  equations  for  failure  pressure  again.  Usually  only  two  or  three  iterations  are  needed  for 
satisfactory  convergence  of  the  failure  pressure  [10]. 

4.2  Asymmetric  Buckling 

Asymmetric  buckling,  or  lobar  buckling,  is  the  collapse  of  the  shell  between  adjacent 
rings  characterized  by  circumferential  lobes  extending  partially  around  the  periphery  of  the 
cylinder.  As  discussed  previously,  this  failure  mode  normally  occurs  when  the  slenderness  ratio 
(A)  is  greater  than  1.14.  Asymmetric  buckling  can  also  occur  when  the  cylinder  shell  is 

relatively  thin  and  the  ring-stiffeners  are  widely  spaced.  In  1929,  von  Mises  first  proposed  a 
solution  to  the  buckling  of  non-reinforced  cylinders  under  hydrostatic  pressure.  He  assumed 
sinusoidal  displacements  in  the  axial  and  circumferential  directions  to  allow  solving  a  set  of 


linearized  partial  differential  equations.  The  equations  represented  the  elastic  action  of  the  shell 
[1].  von  Mises  eventually  obtained  the  following  well  know  equation  for  the  buckling  pressure: 


Where  L  is  the  unsupported  shell  length  between  ring-frames  (L  =  Lf  -  b).  In  this  equation,  the 
buckling  pressure  is  dependent  on  the  number  of  circumferential  lobes  (n),  which  is  an  integer 
value.  To  arrive  at  the  correct  failure  pressure,  an  iterative  process  is  required  varying  n  until  the 
lowest  pressure  is  determined. 

Another  approach  to  minimizing  the  failure  pressure  in  equation  (10)  is  to  solve  it 
analytically,  and  thus  find  an  expression  for  failure  pressure  that  is  independent  of  n.  In  1933, 
Widenburg  solved  this  equation  that  resulted  in  the  Widenburg  approximation  shown  below  [1]: 


P  cLB  := 


2.42  E- 
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Test  data  shows  that  buckling  pressures  determined  by  the  use  of  equation  (11)  differ  by  no  more 
than  about  3.5  percent  from  those  found  from  equation  (10)  [1].  Because  of  its  ease  of 
calculation  and  good  results,  the  Widenburg  approximation  (11)  is  generally  accepted  in  the 
reference  material  as  the  best  method  to  calculate  asymmetric  buckling,  and  was  therefore  used 
for  this  analysis. 


4.3  General  Instability 

General  instability  is  characterized  by  the  failure  of  both  the  shell  and  the  ring-stiffeners. 
A  cylinder  normally  fails  by  general  instability  when  the  rings  are  relatively  “light”  or  “weak”  in 
comparison  to  the  shell,  and  the  cylinder  is  long  [1].  General  instability  can  initiate  in  either  the 
elastic  or  inelastic  stress  region,  but  die  final  configuration  is  in  the  plastic  range  of  the  material. 
Elastic  general  instability  is  die  mode  covered  by  the  available  literature  and  is  addressed  in  this 
thesis.  Inelastic  general  instability  has  been  studied  mainly  by  government  laboratories  and 
organizations.  Most  of  the  material  is  classified  in  nature  and  therefore  not  covered  in  this 
analysis. 

The  first  analysis  of  general  instability  was  conducted  by  Tokugawa  in  1929  [1].  His 
methodology  was  based  on  the  method  of  “split  rigidities”,  where  he  considered  die  failure  of  the 
ring  and  shell  separately  and  summed  the  combined  pressures  [2].  In  the  1940’s  Kendrick  used  a 
strain  energy  method,  with  good  results,  to  determine  the  failure  pressure.  Kendrick’s  solution 
was  rather  complicated  though,  and  in  1954  Bryant  used  a  simpler  strain  energy  method  and 
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developed  an  equation  that  produced  nearly  the  same  results  [2].  The  Bryant  equation  is  used  in 
this  analysis  and  is  shown  below: 


(12) 


Where: 


In  the  Bryant  equation  (12),  the  first  term  corresponds  to  the  shell  failure  and  the  second 
term  to  the  ring  failure,  similar  to  the  “split  rigidity”  used  by  Tokugawa  [2].  The  moment  of 
inertia  (Ie)  used  is  that  of  the  combined  section  of  one  ring  plus  an  “effective”  length  (Le)  of  the 
adjacent  shell.  This  effective  length  term  has  received  significant  attention  over  the  years.  For 
the  purpose  of  this  analysis,  Le  was  calculated  using  the  equation  from  Pulos  and  Salerno  shown 
below  [10]: 


cosh(8)  -  cos  (e)^ 
sinh(G)  +  sin(0)  ) 


(13) 


Where: 


In  order  to  determine  the  failure  pressure  for  general  instability  (equation  12),  the  number 
of  circumferential  lobes  (n)  must  be  varied  to  find  the  number  that  minimizes  the  failure 
pressure. 

The  calculations  and  results  of  the  analytic  solutions  are  provided  in  Appendix  A. 


CHAPTER  5:  Numerical  Solution 

Numerical  analysis  methods  are  very  widely  used  in  engineering  design,  and  are 
employed  extensively  in  the  analysis  of  solids,  structures,  and  fluids.  With  the  advent  of  die 
digital  computer,  the  effectiveness  and  general  applicability  of  this  form  of  engineering  analysis 
was  finally  made  practical.  The  tools  available  for  numerical  analyses  cover  a  wide  range  of 
applicability  from  providing  good  first-order  design  predictions  to  detailed  stress  analyses.  With 
the  increasing  fidelity  of  the  numerical  analysis  tool,  the  cost  of  use  (time  and  money)  also 
increases.  Some  numerical  analysis  tools  are  good  for  initial  design  estimates  and  predictions, 
allowing  the  designer  to  easily  input  data  and  test  several  model  variations.  Other  numerical 
analysis  tools  involve  finite  element  analysis  and  provide  detailed  local  stress  evaluation  of 
structures,  but  involve  complicated  models  that  are  time  consuming  to  develop  and  analyze.  The 
submarine  designer  must  consider  the  benefits  and  applicability  of  the  various  numerical  analysis 
tools  and  determine  which  one  is  appropriate  for  the  particular  stage  of  the  design  process. 

For  this  thesis,  a  numerical  analysis  tool  was  used  for  comparison  to  the  analytic  and 
classification  society  solutions  for  the  failure  pressure  and  failure  mode  of  the  test  cylinders.  As 
a  result,  the  numerical  analysis  tool  was  selected  based  on  its  ease  of  use  and  applicability  for 
cylindrical  structures.  The  tool  was  intended  to  be  used  for  initial  design  predictions  and  not 
local  stress  analysis.  This  thesis  was  not  intended  to  make  comparisons  of  different  numerical 
analysis  tools,  rattier  to  select  one  tool  and  compare  it  to  other  solutions  using  different 
methodology. 


5.1  MAESTRO™  Overview 


The  Method  for  Analysis,  Evaluation,  and  Structural  Optimization  (MAESTRO™)  is  a 
finite  element-based  computer  analysis  tool,  designed  specifically  to  facilitate  the  modeling  of 
ocean  engineering  structures,  including  ships  and  ring-reinforced  cylinders.  MAESTRO™  was 
selected  as  the  tool  for  determining  the  numerical  solution  and  has  the  following  features: 

1)  Rationally-based  analysis  tool,  in  that  it  is  based  on  the  limit-state  approach  to  structural 
design  as  described  in  reference  [7]. 

2)  Capable  of  modeling  virtually  an  entire  structure;  for  a  pressure  hull  of  a  submarine,  this 
includes  the  hull  plating,  frames,  kingframes,  and  bulkheads  to  almost  any  level  of  detail. 

3)  Capable  of  modeling  virtually  any  load  or  combination  of  loads. 

4)  Can  be  operated  in  analysis,  evaluation,  or  optimization  modes. 

The  program’s  underlying  theory  and  detailed  description  of  its  principal  features  are  given  in 
reference  [7],  which  constitutes  the  Theoretical  Manual  for  the  program. 

The  basic  units  of  structural  modeling  are  principal  ship  structural  members  such  as 
beams,  stiffened  panels,  or  girders.  In  order  to  have  an  efficient  interaction  for  the  finite  element 
analysis,  the  elements  used  by  MAESTRO™  are  in  most  cases  the  same  as  the  principal  ship 
structural  members  [1 1].  Elements  are  combined  to  make  strakes  that  are  further  grouped  into 
modules.  A  module  is  a  portion  of  the  structure  being  modeled  that  has  regularly  spaced  sections 
and  local  element  dimensions  that  are  similar;  that  is,  plate  thickness  and  flange  and  web  widths 
and  thicknesses.  Modules  are  then  combined  together  to  create  the  complete  mathematical 
model.  The  mathematical  model  is  meshed  using  several  finite  element  types  discussed  in  detail 
in  reference  [1 1].  MAESTRO™  uses  an  interactive  graphics  program,  MAESTRO™  Modeler, 


to  facilitate  the  creation  of  the  structural  model  and  the  input  file  for  analysis  by  the  source 
program. 

5.2  Cylindrical  Models 

MAESTRO™  is  particularly  useful  for  the  submarine  designer  due  to  its  ability  to 
analyze  cylindrical  structures.  Strakes  can  be  identified  as  part  of  a  complete  cylinder  (360 
degrees)  with  the  curvature  (segment  height,  H)  being  defined  at  the  strake  level  by  the  following 
equation: 

H  :=  R.f  1  -  cosf — ^ 

l  \2))  (14) 

Where: 

0  is  the  strake’s  sector  angle 

When  the  cylinder  option  is  used,  it  implies  that  the  module  includes  one  complete  cylinder  (or 
half  cylinder)  and  that  all  of  the  strakes  are  part  of  that  cylinder.  For  strakes  identified  in  this 
manner,  calculations  are  made  to  determine  the  proximity  to  failure  modes  similar  to  those 
defined  in  Chapter  1.1 

For  this  thesis,  cylindrical  models  were  developed  for  the  test  cylinders  with  known 
dimensions  and  failure  modes  and  pressures.  The  MAESTRO™  Modeler  was  used  to 
graphically  create  the  models  (and  input  file),  while  the  MAESTRO™  (version  8.5)  solver  was 
used  to  perform  the  numerical  analysis. 

5.3  Failure  Modes  Evaluated 

MAESTRO™  uses  limit  states  (or  adequacy  parameters)  for  determining  proximity  to 
failure  for  various  structural  members.  When  using  the  cylinder  feature  in  MAESTRO™, 


specific  calculations  are  invoked  which  replace  some  of  the  limit  state  analyses  with  three  types 
of  cylinder  collapse:  bay  buckling,  general  buckling,  and  local  buckling.  Calculations  for  these 
cylinder  failure  modes  are  based  on  API  (Bull  2U,  1987  edition).  A  detailed  discussion  of  API 
(Bull  2U)  is  provided  Chapter  6. 

The  failure  mode  and  pressure  for  each  test  cylinder  was  determined  by  varying  the 
submergence  pressure  (load)  applied  to  each  model.  Once  one  of  the  limit  states  was  exceeded, 
the  pressure  was  recorded  as  the  failure  pressure  along  with  the  corresponding  failure  mode. 
Results  of  the  numerical  analysis  conducted  using  MAESTRO™  are  provided  in  Appendix  B. 


CHAPTER  6:  Classification  Society  Design  Rules 

The  two  classification  society  design  rules  that  were  utilized  were  the  American 
Petroleum  Institute  and  Det  Norske  Veritas.  These  mles  were  chosen  for  their  availability,  their 
widespread  use  around  the  world,  and  their  coverage  of  the  specific  geometries  of  the 
experimental  cylinders.  Additionally,  API  and  DNV  use  semi-empirical  formulations  that  could 
be  contrasted  to  design  rules  that  use  strictly  closed-form  analytical  equations. 

6.1  American  Petroleum  Institute  (API) 

The  API  design  rules,  as  delineated  in  the  Bulletin  on  Stability  Design  of  Cylindrical 
Shells  [5],  gives  a  brief  and  conservative  approach  for  determining  the  failure  pressures  and 
stresses  for  each  of  the  failure  modes  considered.  Since  the  API  design  rules  are  used  for  many 
different  types  of  marine  structures,  it  accounts  for  several  different  stiffener  and  stringer 
geometries.  The  appropriate  geometry  for  use  in  submarine  design  is  classified  as  a  “ring- 
stiffener”  geometry.  Under  the  ring-stiffener  geometry,  API  (Bull  2U)  addresses  the  following 
buckling  modes:  Local  Shell  Buckling,  General  Instability,  Local  Stiffener  Buckling,  and 
Column  Buckling.  For  comparison  purposes,  local  shell  buckling  and  general  instability 
described  in  [5]  are  the  same  as  asymmetric  buckling  and  general  instability  described  in  Chapter 
1.1,  respectively.  Column  buckling  is  of  concern  for  large  risers  used  to  support  axial  loads 
while  local  stiffener  buckling  is  of  concern  for  designs  with  “light”  stiffeners.  Because  column 
buckling  and  local  stiffener  buckling  are  not  of  concern  for  cylinders  of  the  overall  size  and 
dimensions  used  in  this  analysis,  these  buckling  modes  were  not  considered. 

The  buckling  strength  formulations  presented  in  this  bulletin  are  based  upon  classical 
linear  theory  that  is  modified  by  reduction  factors  to  account  for  the  effects  of  imperfections, 


boundary  conditions,  nonlinear  material  properties  and  residual  stresses.  The  reduction  factors 
are  determined  from  empirical  data  on  shells  of  representative  size  and  initial  imperfections  [5]. 

6.1.1  Limitations  and  Applicability 

API  (Bull  2U)  contains  semi-empirical  formulations  for  evaluating  the  buckling  strength 
of  stiffened  cylindrical  shells.  The  empirical  data  for  these  formulations  was  obtained  through 
numerous  tests  of  ring-stiffened  cylindrical  shells.  As  a  result,  the  failure  modes  and  pressures 
predicted  are  very  dependent  upon  having  cylinders  similar  in  size  to  those  used  for  the  empirical 
data. 

API  (Bull  2U)  is  applicable  to  shells  that  are  fabricated  from  steel  plates  where  the  plates 
are  cold  or  hot  formed  and  joined  by  welding,  and  stiffeners  are  to  be  uniformly  spaced.  It  is 
intended  for  design  and  review  of  large  diameter  cylindrical  shells,  typically  identified  as  those 
with  diameter  to  shell  thickness  (D/t)  ratios  greater  than  300  but  less  than  2000.  A  minimum 
shell  thickness  of  3/16  inches  is  allowed  with  a  limit  of  shell  radius  to  shell  thickness  (R/t)  ratio 
of  less  than  or  equal  to  1000.  Most  of  the  material  used  for  empirical  tests  had  yield  strengths 
between  36  ksi  and  100  ksi  [5]. 

6.1.2  Local  Shell  Buckling 

The  failure  pressure  for  local  buckling  mode  (asymmetric  buckling)  was  determined  by 
first  solving  for  the  theoretical  failure  pressure  for  local  buckling  (peL>  of  the  cylindrical  shell. 
Once  the  theoretical  failure  pressure  was  known,  reduction  factors  were  applied  that  account  for 
fabrication  tolerances  (o^l)  and  plasticity  reduction  (q)  for  nonstress  relieved  shells  [5].  The 
equation  for  determining  the  failure  pressure  for  local  buckling  mode  is  shown  below: 


The  theoretical  failure  pressure  formulation  was  based  on  an  equation  that  was  derived  from  von 
Mises  equation  (10).  The  theoretical  failure  pressure  is  a  smooth  lower  bound  curve  of  (10) 
which  was  obtained  by  letting  the  number  of  circumferential  lobes  (n)  be  a  noninteger  [5].  The 
equations  for  the  theoretical  failure  pressure  are  shown  below: 


Where: 


P  eL := 


5.08 


a  118  n  c 

A  m  +  °-5 


—•&(—) 

Am 


if  M  x  >  1.5  a  A  m  <  2.5 


if  2.5  <  Am<0.104-^yJ 


6.688-Cp-  1061'e(^j  j 


if  0.208  <Cp  <2.85 


if 


Cp  >2.85 


(16) 


Mx" 


■JRA 


Am:=Mx-  1.17+  1.068k 


Cp: 


(2Aj 

(?) 


k  =  0.5  for  external  hydrostatic  pressure 

Imperfection  factors  (ool)  are  generally  assigned  a  constant  value  of  0.8  for  fabrication  processes 
that  meet  the  specifications  given  in  [5].  The  plasticity  reduction  factors  are  applied  using  the 
following  equations: 


11  := 


Where: 


1  if  Ac  <  0.55 
0.45  'N 

- +  0.18  if  0.55  <Ac<,  1.6 

VAc  J 


- : -  if  1.6  <  Ac  <  6.25 

1  +  1.15-Ac 

—  if  Ac  S:  6.25 
Ac 


Fy=  yield  strength  of  the  material 
FreLr  =  elastic  buckling  stress 


(17) 


6.1 3  General  Instability 


The  failure  pressure  for  general  instability  mode  was  calculated  by  first  determining  the 
theoretical  elastic  failure  pressure  for  general  instability  (peo)  of  the  ring-stiffened  shell.  Once 
the  elastic  failure  pressure  was  determined,  reduction  factors  were  applied  that  account  for 
fabrication  tolerances  (dec)  and  plasticity  reduction  (x\)  for  nonstress  relieved  shells,  in  a  manner 


similar  to  Chapter  6.1.2.  The  equation  for  determining  the  failure  pressure  for  general  instability 


mode  is  shown  below: 

PcGr:-1  °eaPec{") 


(18) 


The  theoretical  failure  pressure  formulation  was  based  on  the  Bryant  equation  (12),  where  the 
failure  pressure  was  a  function  of  the  number  of  circumferential  lobes  (n),  which  must  be  varied 
to  determine  the  minimum  pressure  value.  The  equation  for  the  theoretical  failure  pressure  is 
shown  below: 
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(19) 


6.2  Det  Norske  Veritas  (DNV) 

The  DNV  design  rules,  as  delineated  in  the  Buckling  Strength  of  Shells,  Recommended 
Practice  DNV-RP-C202  [6],  treats  the  buckling  stability  of  shell  structures  based  On  the  load  and 
resistance  factor  design  format  (LRFD).  The  methods  used  in  [6]  are  considered  semi-empirical. 
The  reason  for  basing  the  design  on  semi-empirical  methods  is  that  the  agreement  between 
theoretical  and  experimental  buckling  loads  for  some  cases  has  been  found  to  be  non-existent. 
This  discrepancy  is  due  to  the  effect  of  the  geometric  imperfections  and  residual  stresses  in 
fabricated  structures.  Actual  geometric  imperfections  and  residual  stresses  do  not  in  general 
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appear  as  explicit  parameters  in  the  expressions  for  buckling  resistance.  This  means  that  the 
methods  for  bucking  analysis  are  based  on  an  assumed  level  of  imperfections.  For  DNV,  this 
tolerance  level  is  specified  in  DNV  OS-C401;  Fabrication  and  Testing  of  Offshore  Structures. 
Since  the  DNV  design  rules  are  used  for  many  different  types  of  marine  structures,  they  account 
for  several  different  stiffener  and  stringer  geometries.  The  appropriate  geometry  for  use  in 
submarine  design  is  classified  as  a  “ring-stiffened”  geometry.  Under  the  ring-stiffened 
geometry,  DNV  (RP-C202)  addresses  the  following  buckling  modes:  Shell  Buckling,  Panel 
Ring  Buckling,  and  Column  Buckling.  For  comparison  purposes,  shell  buckling  and  panel  ring 
buckling  described  in  [6]  are  the  same  as  asymmetric  buckling  and  general  instability  described 
in  Chapter  1.1,  respectively.  Column  buckling  is  of  concern  for  large  risers  used  to  support  axial 
loads.  Because  column  buckling  is  not  of  concern  for  cylinders  of  the  overall  size  and 
dimensions  used  in  this  analysis,  this  buckling  mode  was  not  considered. 

6.2.1  Limitations  and  Applicability 

Similar  to  API  (Bull  2U),  DNV  (RP-C202)  contains  semi-empirical  formulations  for 
evaluating  the  buckling  strength  of  stiffened  cylindrical  shells.  In  the  case  of  the  DNV  design 
rules  though,  no  specific  limitations  were  placed  on  the  size  or  dimensions  of  the  cylinders.  The 
only  specified  requirement,  assumes  the  edges  of  the  cylinder  are  effectively  supported  by  ring 
frames,  bulkheads  or  end  closures  [6].  Neither  empirical  data  nor  experimental  results,  used  to 
derive  the  equations  for  cylinder  buckling  were  provided  in  [6]. 

6.2.2  Shell  Buckling 

The  failure  pressure  for  shell  buckling  (asymmetric  buckling)  was  determined  by  first 
calculating  the  characteristic  buckling  strength  of  the  shell  (fks).  The  characteristic  buckling 


strength  was  then  divided  by  a  material  factor  (yM)  to  determine  the  design  shell  buckling 
strength  (fksd)  shown  below: 

^ks 


Where: 


[ksd  • 


m 


(20) 


^  ks  ‘ 


Ym:= 


i 


1  +  X 


ssquare 


1.15  if  ?ls<0.5 


(0.85+0.6Xs)  if  0.5<XS  <,  1.0 


1.45  if  >1.0 


To  solve  (20),  the  reduced  shell  slenderness  (Xs)  must  first  be  defined.  Xg  is  a  combination  of  the 
shell  stresses  and  elastic  buckling  strength.  For  the  current  analysis,  the  design  stresses 
associated  with  bending  and  shear  were  neglected  since  they  were  not  present  in  the  test 
cylinders  analyzed.  The  equation  for  reduced  shell  slenderness  is  shown  below: 

(21) 


F  f 
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'  s_square 
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j_sd 


sd  ^Lh  sd 


v  fEa 
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Where: 


CTj_sd :-^CTa_sd  CTa_sd‘CTh_sd  +  CTh_sd  ) 
-press  R 


a  sd  • 


2t 


design  equivalent  von  Mises  stress 


design  axial  stress 


rh  sd' 


-prcss-R 


a  + 1 


design  circumferential  stress 
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CSBn  E  f  t  ^ 
fc  := — 7 - r-  —  1 


:.(l-v2)lLrJ 


elastic  buckling  strength  of  shell 


c®:=4+(v) 

The  remaining  constants  and  definitions  are  provided  in  [6]. 


reduced  buckling  coefficient 


Using  the  above  equations,  the  design  shell  buckling  strength  (fksd)  can  be  determined  for 
a  specific  submergence  pressure.  The  shell  buckling  strength  was  then  compared  to  the  design 
equivalent  von  Mises  stress  (aj  sd).  If  Oj_Sd<fksd  then  the  cylinder  should  not  fail  by  shell 
buckling.  To  determine  the  pressure  at  which  the  shell  would  fail,  oj  sd>fksd,  an  iterative  process 
was  used.  Submergence  pressure  was  increased  in  step  intervals  while  recalculating  (20)  and 
(21).  Once  the  limit  condition  was  exceeded,  a  failure  pressure  for  shell  buckling  was 
determined. 

6.2.3  Panel  Ring  Buckling 

Failure  by  panel  ring  buckling  (general  instability)  was  determined  by  evaluating  both  the 
cross  sectional  area  of  a  ring  frame  and  the  effective  moment  of  inertia  of  a  ring  frame.  To 
ensure  the  ring  frame  would  not  fail  prematurely  the  cross  sectional  area  of  a  ring  frame 
(exclusive  of  the  effective  shell  plate  flange)  should  not  be  less  than  AReq,  defined  by: 


A  Req  ^ 


—  +  0.06^L,t 


Where: 


Zi  := — Jl  -v 
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If  a  ring-stiffened  cylinder,  or  a  part  of  a  ring-stiffened  cylinder,  is  effectively  supported 
at  the  ends,  a  refined  calculation  of  moment  of  inertia  (Ih)  is  used  by  DNV  (RP-C202)  for 
calculating  the  capacity  of  the  ring  frame.  Using  an  initial  geometry,  an  effective  moment  of 
inertia  of  die  combined  ring  frame  and  shell  (Ie)  can  be  calculated.  The  value  for  Ie  is  also 
implicit  in  the  procedure  for  calculating  the  buckling  capacity  of  the  panel  and  ring. 

When  a  ring-stiffened  cylinder  is  subjected  to  external  pressure,  the  ring-stiffeners  should 

satisfy: 


1  + 


P  sd  GI^0-75 - t  rf 
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1  characteristic  buckling  strength 

fr  =  characteristic  material  yield  strength  (Fy) 


reduced  column  slenderness 
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effective  radius  of  gyration  or  ring  frame 

equivalent  length  of  shell  plating 
initial  out-of-roundness  parameter 
effective  moment  of  inertia 

curvature  parameter 
Ar  =  cross-sectional  area  of  ring  frame 
rf  =  radius  of  shell  measured  to  ring  flange 

Using  the  above  equations,  the  maximum  design  external  pressure  can  be  determined. 
Given  known  cylinder  and  ring  frame  dimensions,  values  can  be  substituted  in  equation  (23)  and 
the  associated  calculations  to  determine  a  maximum  design  external  pressure.  The  maximum 
design  pressure  is  then  considered  the  failure  pressure  for  the  panel  ring  failure  (general 
instability)  mode. 

The  results  and  calculations  of  the  DNV  (RP-C202)  analysis  are  provided  in  Appendix  D. 
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CHAPTER  7:  Results 


The  analytic  solution,  numerical  solution,  and  classification  society  design  solution  were 
all  compared  against  test  data  collected  from  experiments  conducted  by  the  United  States  Navy 
in  support  of  submarine  design.  Each  solution  method  was  used  to  determine  a  failure  pressure 
and  failure  mode  for  each  of  the  test  cylinders.  The  resulting  predictions  were  then  compared  to 
the  experimental  results.  Of  primary  interest  was  the  agreement  between  the  predicted  and 
experimental  mode  of  failure,  followed  by  the  accuracy  of  the  predicted  failure  pressure  when 
compared  to  the  actual  failure  pressure. 

7.1  NAVSEA  Test  Cylinders 

The  test  data  was  provided  by  the  Naval  Sea  Systems  Command  Submarine  Structural 
Integrity  Division  (NAVSEA  05P2).  Data  was  provided  for  four  test  cylinders  that  were  selected 
to  cover  the  range  of  examined  failure  modes.  The  cylinder  diameter  to  thickness  ratios  (D/t) 
were  from  1 12  to  198.  Two  of  the  cylinders  had  internal  stiffeners  while  the  other  two  cylinders 
had  external  stiffeners.  All  four  test  cylinders  had  built-up  end  stiffeners  with  a  combination  of 
narrower  spacing  and  /  or  larger  stiffener  dimensions  than  the  uniform  section  of  stiffeners.  The 
end  stiffeners  were  designed  to  prevent  shell  yielding  in  the  end  bays  due  to  increased  stress 
levels  associated  with  the  boundary  conditions.  It  was  estimated  that  without  the  end  stiffeners  a 
4-5%  reduction  in  axisymmetric  yielding  pressure  could  occur  [12].  Neither  the  analytic  solution 
nor  the  classification  society  design  rules  allowed  for  variable  spaced  stiffeners,  therefore  the 
non-uniformities  were  disregarded  and  the  end  stiffeners  were  treated  as  uniform  section 
stiffeners.  The  four  test  cylinders  are  described  below. 


7.1.1  Cylinder  l.d 


Cylinder  1  .d  was  a  machined  cylindrical  shell  with  external  rectangular  stiffeners.  The  • 

material  used  was  high  strength  steel  with  yield  strength  of  80,000  psi.  Figure  5  shows  the 
structure  and  principal  dimensions  of  the  cylinder.  The  boundary  conditions  consisted  of  one 
end  being  fully  fixed  with  the  other  end  having  all  freedoms  fixed  except  for  axial  displacement. 

External  hydrostatic  pressure  was  applied  including  axial  line  load  to  simulate  load  on  the  end 

plate.  The  experiment  tested  the  ability  of  the  analysis  method  to  predict  elastic  shell  bucking  0 

(asymmetric  buckling).  The  experimentally  determined  collapse  pressure  was  633  psi  with 

failure  by  asymmetric  (Lobar)  buckling. 
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7.1.2  Cylinder  l.f 


•  Cylinder  l.f  was  a  cylindrical  shell  with  internal  tee  stiffeners  of  welded  construction. 

The  material  used  was  high  strength  steel  with  yield  strength  of 98,500  psi.  The  boundary 
conditions  consisted  of  4.0  inch  steel  plates  attached  with  full  fixity  to  the  end  of  the  adaptor  ring 
on  the  model.  External  hydrostatic  pressure  was  applied.  This  test  cylinder  was  used  to  predict 
failure  by  elastic  general  instability.  There  was  no  experimental  elastic  collapse  pressure; 

^  therefore  the  critical  pressure  was  calculated  by  two  separate,  reliable  analysis  programs  with  the 

results  being  4858  psi  (with  3  waves)  and  4953  psi  (with  3  waves).  The  test  cylinder  actually 
failed  by  inelastic  general  instability  at  a  pressure  of  2200  psi.  Figure  6  shows  the  structure  and 
®  dimensions  of  the  test  cylinder. 

Figured:  Test  Cylinder  l.f  Structural  Dimensions 
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7.1.3  Cylinder  2.a 


Cylinder  2.a  was  a  machined  cylindrical  shell  with  external  tee  stiffeners.  The  material 
used  was  high  strength  steel  with  yield  strength  65,500  psi.  Figure  7  shows  the  structure  and 
dimensions  of  the  test  cylinder.  The  boundary  conditions  consisted  of  end  closures  made  of  3.0 
inch  steel  plates  attached  to  the  idealized  adaptor  ring  with  full  fixity.  External  uniform 
hydrostatic  pressure  was  applied  to  the  model.  The  cylinder  tests  die  ability  of  the  analysis 
methods  to  predict  end  bay  failure  (shell  collapse  influenced  by  end  bay  design).  Specifically 
this  model  provides  an  example  of  axisymmetric  buckling.  The  experimental  collapse  pressure 
was  determined  to  be  921  psi  by  axisymmetric  collapse  in  the  second  bay  from  the  adaptor  ring. 

Figure  7:  Test  Cylinder  2.a  Structural  Dimensions 


Typical  Frames 


End  Frames 


7.1.4  Cylinder  2.c 


Cylinder  2.c  was  a  fabricated  cylinder  with  internal  ring-stiffeners.  The  base  material 
used  was  high  strength  steel  with  yield  strength  of  157,000  psi.  Figure  8  shows  the  structure  and 
dimensions  of  the  cylinder.  The  shell  was  cold  rolled  and  fabricated  with  a  deliberate  out-of- 
roundness  imperfection.  The  frames  were  built-up.  The  frame  web  material  was  base  metal,  and 
die  frame  flanges  were  cold  rolled.  The  boundary  conditions  consisted  of  one  end  being  fully 
fixed  with  the  other  end  having  all  freedoms  except  axial  displacement.  External  uniform 
hydrostatic  pressure  with  an  axial  end  load  to  simulate  end  plate  loading  was  applied.  The  test 
cylinder  was  used  to  predict  the  inelastic  general  instability  failure  mode  and  to  model  out-of- 
roundness  imperfections.  In  the  current  analysis  the  out-of-roundness  was  not  considered,  and 
inaccuracies  in  predicted  results  were  expected.  The  collapse  pressure  was  experimentally  found 
to  be  3640  psi  in  two  circumferential  waves  in  an  inelastic  general  instability  mode. 


Figure  8:  Test  Cylinder  2.c  Structural  Dimensions 
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7.2  Calculation  to  Experimental  Comparison 


The  dimensions  of  each  of  the  test  cylinders  were  used  to  calculate  a  predicted  failure 
mode  and  failure  pressure  for  the  analytic  solution,  numerical  solution,  and  classification  society 
design  rules  (Appendix  A-D).  Table  1  shows  the  comparison  of  die  predicted  solutions  to  the 
experimental  results.  The  table  displays  the  calculated  failure  mode  and  failure  pressure  for  each 
test  cylinder. 

Table  1:  Comparison  of  Predicted  Failure  Mode  and  Pressure  to  Experimental  Data 


HI 

l.d 

r  CyTTf  “ 

Cyl 

2.a 

Cvl 

2.c 

Mode 

Pressure 

Mode 

Pressure 

Mode 

Pressure 

Mode 

Experiment 

633 

L 

2200 

Gl 

921 

AX 

3640 

Gl 

623 

L 

2176 

AX 

885 

AX 

4137 

AX 

API  (Bull  2U) 

447 

L 

1838 

Gl 

710 

L 

3784 

Gl 

HUH 

164 

L 

1089 

Gl 

487 

Gl 

2457 

Gl 

MAESTRO 

567 

L 

1920 

Gl 

797 

L 

4167 

Gl 

Key:  L  Asymmetric  (Lobar)  Buckling 

AX  Axisymmetric  Yielding 

Gl  General  Instability 


API,  DNV,  and  MAESTRO  only  address  L  and  Gl  for  ring-stiffened  cylinders 


The  comparative  analysis  for  each  of  the  test  cylinders  is  discussed  in  the  following  sections. 

7.2.1  Cylinder  l.d  Results 

Cylinder  1  .d  was  designed  to  test  the  ability  of  design  tools  to  predict  asymmetric 
buckling  failure.  There  was  excellent  agreement  between  the  predicted  failure  modes  and  the 
experimentally  determined  failure  mode.  Asymmetric  buckling  was  expected  since  the 


slenderness  ratio  was  201,  much  greater  than  the  breakpoint  (1.14)  between  asymmetric  and 
axisymmetric  failure.  Specific  comparisons  for  cylinder  l.d  are  shown  in  Table  2  and  Figure  9. 


Table  2:  Predicted  Failure  Pressures  and  Modes  for  Cylinder  l.d 
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Figure  9:  Comparison  of  Predicted  Failure  Pressures  for  Cylinder  l.d 

The  predicted  failure  pressures  covered  a  large  range  of  values.  The  closed-form  analytic 
solution  was  very  close  to  the  experimental  value,  being  only  1.6%  below  the  critical  pressure. 
The  classification  society  design  rules  significantly  under  predicted  the  failure  pressure.  In  the 


47 


case  of  API  (Bull  2U),  the  discrepancy  is  attributable  to  the  reduction  factors  applied  for  cylinder 
imperfections  and  residual  stresses  due  to  fabrication.  Test  cylinder  1  .d  was  machined,  therefore  • 

these  reduction  factors  would  not  be  applicable.  DNV  (RP-C202)  is  based  on  a  semi-empirical 
solution  and  the  results  are  closely  tied  to  the  specific  cylinders  used  to  develop  the  empirical 
relationships.  Cylinder  1  .d  is  relatively  small  and  therefore  is  not  modeled  well  by  the  DNV  * 

design  rules.  The  MAESTRO™  solution  under  predicted  the  failure  pressure  by  approximately 

10%. 


7.2.2  Cylinder  l.f  Results 

Cylinder  1  .f  was  fabricated  to  test  the  ability  of  design  tools  to  predict  general  instability.  ^ 

The  cylinder  failed  under  experiment  at  2200  psi  by  general  instability.  There  was  excellent 
agreement  between  the  predicted  failure  modes  and  the  experimentally  determined  failure  mode. 

Specific  comparisons  for  cylinder  1  .f  are  shown  in  Table  3  and  Figure  10.  • 

Table  3:  Predicted  Failure  Pressures  and  Modes  for  Cylinder  l.f 
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Cylinder  l.f  Comparisons 


Design  Tool 


Figure  10:  Comparison  of  Predicted  Failure  Pressures  for  Cylinder  l.f 

All  solutions  except  the  closed-form  analytic  solution  predicted  that  the  cylinder  would 

fail  by  general  instability.  The  analytic  solution  predicted  that  the  failure  would  be  by 

axisymmetric  yielding  at  a  pressure  that  was  very  close  to  the  experimental  failure  pressure  (- 

* 

1.1%).  In  the  case  of  API  (Bull  2U),  a  form  of  the  Bryant  equation  was  used  to  calculate  the 
failure  pressure  for  general  instability,  and  the  difference  from  experimental  value  is  attributed  to 
the  difference  between  reduction  factors  and  actual  imperfections  and  residual  stresses.  DNV 
(RP-C202)  again  significantly  under  predicted  the  failure  pressure.  Cylinder  l.f  was  slightly 
larger  than  cylinder  l.d,  but  still  not  on  the  size  scale  of  typical  marine  structures,  contributing  to 
the  difference  found  in  the  DNV  prediction.  The  MAESTRO™  solution  under  predicted  the 
failure  pressure  by  approximately  13%. 


7.2.3  Cylinder  2.a  Results 


Cylinder  2.a  was  fabricated  to  test  the  ability  of  the  design  tools  to  predict  axisymmetric 
yielding.  Of  the  solutions  used  for  the  current  analysis,  only  the  analytic  solution  addresses 
axisymmetric  yielding  and  provides  a  failure  pressure.  Test  cylinder  2.a  was  experimentally 
determined  to  fail  by  axisymmetric  yield  at  921  psi.  Specific  comparisons  for  cylinder  2.a  are 
shown  in  Table  4  and  Figure  11. 

Table  4:  Predicted  Failure  Pressures  and  Modes  for  Cylinder  2.a 


failure  Pressure 
(psi) 

Failure 

Mode 

- %Trom - 

Experiment 

Experiment 

921 

AX 

— 

Analytic  Solution 

885 

AX 

-3.9 

API  (Bull  2U) 

710 

L 

-22.9 

DNV  (RP-C202) 

487 

Gl 

-47.1 

MAESTRO 

797 

L 

-13.5 

Figure  11:  Comparison  of  Predicted  Failure  Pressures  for  Cylinder  2.a 


The  predicted  failure  modes  for  cylinder  2.a  varied  significantly  between  the  design  tools 
used.  Axisymmetric  yield  was  expected  since  the  slenderness  ratio  was  0.698,  much  less  than 
the  breakpoint  (1.14)  between  asymmetric  and  axisymmetric  failure.  In  the  case  of  the  analytic 
solution,  the  correct  failure  mode  was  predicted  and  the  failure  pressure  was  within  in  4%  of 
experimental  value,  showing  a  good  correlation.  API  (Bull  2U)  and  DNV  (RP-C202)  do  not 
calculate  an  axisymmetric  yield  failure  pressure.  Most  large  marine  structures  have  a 
slenderness  ratio  greater  than  1.14,  therefore  it  is  presumed  that  API  and  DNV  are  primarily 
concerned  with  these  types  of  structures  and  do  not  consider  axisymmetric  yield  important  for 
their  design  rules.  MAESTRO™  uses  API  (Bull  2U,  1987  edition)  for  its  calculations  and 
therefore  does  not  address  axisymmetric  yield  either. 

7.2.4  Cylinder  2.c  Results 

Cylinder  2.c  was  fabricated  to  test  the  ability  of  the  design  tools  to  predicted  general 
instability  and  model  out-of-roundness  imperfections.  While  none  of  the  design  tools  used 
explicitly  modeled  out-of-roundness  imperfections,  the  classification  society  design  rules  do 
apply  generic  reduction  factors  for  imperfections  and  residual  stresses.  There  was  generally 
good  agreement  between  the  predicted  failure  modes  and  the  experimentally  determined  failure 
mode.  Specific  comparisons  for  cylinder  2.c  are  shown  in  Table  5  and  Figure  12. 

Table  5:  Predicted  Failure  Pressures  and  Modes  for  Cylinder  2.c 


Figure  12:  Comparison  of  Predicted  Failure  Pressures  for  Cylinder  2.c 

All  solutions  except  the  closed-form  analytic  solution  predicted  that  the  cylinder  would 
fail  by  general  instability.  The  analytic  solution  predicted  that  the  failure  would  be  by 
axisymmetric  yield  at  a  pressure  approximately  14%  greater  than  the  actual  failure  pressure.  In 
the  case  of  API  (Bull  2U),  the  predicted  failure  pressure  was  very  close  to  the  actual  failure 
pressure  (4%  over)  with  the  error  attributed  to  the  rough  approximations  used  for  nominal 
imperfection  levels.  DNV  (RP-C202)  once  again  predicted  the  correct  failure  mode,  but  greatly 
under  predicted  the  failure  pressure.  The  MAESTRO™  solution  over  predicted  the  failure 
pressure  by  approximately  15%. 
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7.3  Comparison  to  Previous  Work 

Analysis  was  previously  conducted  by  D.  Price  [4],  comparing  the  failure  modes  and 
pressures  of  the  same  test  cylinders  (cylinders  l.d,  l.f,  2.a,  2.c)  using  different  classification 
society  design  rules  and  guides.  Table  6  shows  the  comparison  of  [4]  to  the  current  analysis. 


Table  6:  Comparison  With  Previous  Work 


Cyl  1.d  ] 

Cyl  l.f  ^ 

Cyl  2.a 

Cyl  2.c  j 

Mode 

Pressure 

Mode 

Pressure 

Mode 

Pressure 

Mode 

Experiment 

633 

■■ 

2200 

Gl 

921 

AX 

3640 

Gl  ! 

Analytic 

623 

H 

2176 

AX 

885 

AX 

4137 

AX 

API  (Bull  2U) 

447 

H 

1838 

Gl 

710 

m 

3784 

Gl 

DNV  (RP-C202) 

164 

H 

1089 

Gi 

487 

Gl 

2457 

Gl 

MAESTRO 

567 

H 

1920 

Gl 

797 

L 

4167 

Gl 

Key: 

L 

Lobar  Buckling 

previous  work 

Note: 


AX  Axisymmetric  Yielding 

Gl  General  Instability 

API,  DNV,  and  MAESTRO  only  address  L  and  Gl  for  ring-stiffened  cylinders 


The  classification  society  design  rules  and  guides  examined  in  [4]  were:  The  American  Bureau 
of  Shipping  (ABS)  ( Rules  for  Building  and  Classing  Underwater  Vehicles,  Systems  and 
Hyperbaric  Facilities,  1990  Edition),  Germanischer  Lloyd  ( Rules  for  Underwater  Technology, 
1998  Edition),  and  Society  of  Naval  Architects  and  Marine  Engineers  (SNAME)  ( Submersible 
Vehicle  Systems  Design,  1990  Edition).  The  design  tools  used  in  [4]  all  predicted  the  same 
failure  mode  as  the  analytic  solution.  This  was  due  primarily  to  the  fact  that  most  of  the 
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solutions  for  these  rules  and  guides  were  based  on  closed-form  solutions  similar  to  those  used  in 
the  analytic  solution.  In  the  case  of  cylinders  l.d  and  2.a,  the  design  rules  and  guides  predicted 
the  correct  failure  mode  for  both  cylinders,  and  failure  pressures  were  all  within  12%  of  the 
experimental  value.  Cylinders  l.f  and  2.c  though,  were  predicted  to  fail  by  axisymmetric  yield 
and  instead  failed  by  general  instability.  This  discrepancy,  when  determining  failure  by  general 
instability,  indicates  there  may  be  some  limitations  with  the  closed-form  analytic  solution 
predicting  a  general  instability  failure. 


CHAPTER  8:  Conclusions 

The  classification  society  design  rules  and  numerical  solution  studied  in  this  thesis  are 
important  tools  for  engineers  and  naval  architects  designing  and  manufacturing  ring-stiffened 
cylindrical  structures  (like  submarines)  subjected  to  external  hydrostatic  pressure.  The  design 
engineer  must  have  confidence  in  his  analysis  tools  and  be  assured  that  they  will  provide 
appropriate,  and  safe  calculations  for  the  structure  under  consideration.  Confidence  can  be 
assured  by  comparison  of  the  calculated  failure  pressure  and  mode  to  those  determined 
experimental  through  physical  model  tests.  This  thesis  attempted  to  provide  that  comparison  by 
using  two  widely  used  classification  society  design  rules  along  with  a  numerical  analysis 
program  to  compare  analytical,  numerical,  and  experimental  results. 

8.1  Comparative  Analysis  Review 

As  discussed  in  Chapter  5,  the  design  rules  and  numerical  methods  had  mixed  results  in 
correctly  predicting  the  failure  mode  and  failure  pressure  of  the  test  cylinders.  In  general,  the 
classification  society  design  rules  and  the  numerical  solution  correctly  predicted  the  failure  mode 
for  the  types  of  failures  addressed  by  these  sources  (asymmetric  buckling  and  general 
instability).  In  the  case  of  axisymmetric  yield  though,  as  found  in  cylinder  2.a,  API  (Bull  2U) 
and  MAESTRO™  predicted  failure  by  asymmetric  buckling  (local  shell  buckling)  and  DNV 
(RP-C202)  predicted  failure  by  general  instability.  Failure  pressures  were  not  predicted  well  by 
the  design  rules  with  API  predicting  failure  at  a  pressure  16%  lower  on  average  than  actual,  and 
DNV  predicting  failure  at  a  pressure  51%  lower  on  average  than  actual.  The  numerical  solution, 
MAESTRO™,  performed  markedly  better  predicting  a  failure  pressure  5.5%  lower  on  average 
than  actual.  In  general,  the  classification  society  design  rules  and  numerical  solution  all  were 


overly  conservative  when  predicting  failure  pressures  for  the  cylinders  fabricated  within  normal 
tolerances.  Cylinder  2.c  was  manufactured  with  a  deliberate  0.105  inch  out-of-roundness  that 
resulted  in  all  the  design  tools,  except  DNV  (RP-C202),  predicting  failure  at  a  pressure  higher 
than  experimental. 

8.2  Agreements  and  Discrepancies 

In  this  thesis  there  were  sixteen  failure  pressures  calculated  (four  failures  for  each  of  four 
test  cylinders).  To  provide  a  useful  comparison  methodology  for  the  design  tools  that  did  not 
address  axisymmetric  yielding,  failure  by  asymmetric  buckling  (lobar)  and  axisymmetric  yield 
were  both  considered  to  be  a  local  failure  of  the  shell  (between  the  stiffeners)  and  considered  a 
similar  failure  mode.  Local  shell  failure  was  then  contrasted  to  general  instability  that  was  a 
failure  of  the  shell  and  ring-stiffener.  Using  these  guidelines  for  comparison,  the  design  tools 
(Analytic,  API,  DNV,  MAESTRO™)  correctly  predicted  81%  of  the  failure  modes.  The  failure 
pressures  calculated  varied  from  experimental  by  1%  to  74%. 

The  closed-form  analytic  solution  correctly  predicted  the  failure  mode  for  two  of  the  four 
test  cylinders  and  predicted  a  failure  pressure  within  4%  for  those  two  cases.  Both  test  cylinders 
that  failed  by  inelastic  general  instability,  were  incorrectly  predicted  to  fail  by  axisymmetric 
yield.  This  indicates  that  the  two  modes  of  failure  are  very  close  together  and  that  the  closed- 
form  solution  has  difficulty  discerning  the  two  failure  modes. 

The  API  (Bull  2U)  solution  correctly  predicted  local  shell  failure  and  general  instability 
in  all  four  of  the  test  cylinders  and  predicted  the  failure  pressure  within  29%  of  the  experimental 
value.  Since  the  API  solution  is  based  on  the  von  Mises  equation  (10)  for  local  shell  buckling 
and  the  Bryant  equation  (12)  for  general  instability,  it  was  expected  that  the  results  would  be 
similar  to  the  analytic  solution.  The  differences  between  the  API  and  experimental  results  can  be 


attributed  to  two  sources.  First,  the  constants  used  in  API  (Bull  2U)  are  derived  from  actual  test 
data.  Being  empirical  in  nature,  these  values  are  closely  correlated  to  the  size  and  dimensions  of 
the  models  tested  to  determine  the  coefficients.  In  the  case  of  the  four  test  cylinders  that  were 
analyzed,  their  dimensions  were  smaller  than  those  used  to  determine  the  empirical  constants. 
Second,  the  reduction  factors  applied  (imperfections  and  residual  stress)  do  not  accurately  model 
the  test  cylinders  (e.g.  two  test  cylinders  were  machined,  where  the  constant  imperfection 
reductions  would  not  be  accurate).  As  a  result,  some  inaccuracies  in  the  predicted  pressures 
were  encountered. 

The  DNV  (RP-C202)  solution  correctly  predicted  local  shell  failure  and  general 
instability  in  three  of  the  four  test  cylinders  but  significantly  under  predicted  the  failure  pressure 
by  33%  to  74%  of  the  experimental  value.  The  differences  between  the  DNV  and  experimental 
results  are  attributed  to  the  semi-empirical  nature  of  the  equations  in  the  rules.  The  constants  and 
equations  used  in  DNV  (RP-C202)  are  empirically  derived  and  therefore  closely  correlated  to  the 
size  and  dimensions  of  the  models  used  to  determine  them.  There  were  no  limitations  on 
applicability  stated  in  the  rules,  but  there  were  large  discrepancies  noted  in  the  predicted  failure 
pressures.  This  appears  to  indicate  the  DNV  rules  are  more  applicable  to  large  marine  structures 
instead  of  smaller  submersible  or  submarine  designs. 

The  MAESTRO™  solution  correctly  predicted  local  shell  failure  and  general  instability 
in  all  four  of  the  test  cylinders  and  predicted  the  failure  pressure  within  15%  of  the  experimental 
value.  The  MAESTRO™  cylinder  analysis  function  is  based  on  the  design  rules  from  API  (Bull 
2U,  1987  edition),  and  was  therefore  expected  to  provide  results  similar  to  those  from  the  design 
rules.  The  MAESTRO™  solution  was  somewhat  more  accurate  than  the  design  rules  because  it 
allows  for  the  use  of  non-uniform  stiffeners  and  spacing.  Additionally,  the  cylinder  can  be 


modeled  in  whole  including  the  end  plates  and  fixturing.  The  added  flexibility  of  MAESTRO™ 
appeared  to  improve  the  accuracy  of  the  predicted  failure  pressures  over  the  other  two 
classification  society  design  rules. 

8.3  Applications  of  the  Results 

The  various  design  rules  and  numerical  methods  studied  are  promulgated  to  ensure  a  safe 
and  adequate  design  of  ring-stiffened  cylinders  for  use  under  external  hydrostatic  pressure.  The 
safe  design  requires  a  high  degree  of  certainty  that  the  cylinder  will  not  fail  under  the  worse  case 
anticipated  conditions.  The  comparison  of  the  closed-form  analytic  solution,  the  classification 
society  design  rules  (API  and  DNV),  and  the  numerical  solution  (MAESTRO™)  to  experimental 
results  allows  a  designer  to  have  a  good  understanding  of  the  strengths  and  limitations  of  each 
analysis  method.  The  ability  to  compare  design  predictions  would  be  useful  in  judging  the  initial 
feasibility  of  a  design  and  also  the  case  where  a  specific  design  is  subject  to  more  than  one 
classification  society. 

The  solutions  and  comparisons  in  this  thesis  should  in  no  way  be  used  as  detailed  design 
tools  for  construction  of  ring-stiffened  cylinders.  Instead,  a  useful  methodology  for  early  stage 
design  can  be  developed  from  these  comparisons.  First,  an  initial  design  can  be  quickly 
developed  using  closed-form  analytic  solutions  similar  to  those  in  Appendix  A.  Next,  the  size 
and  dimensions  of  the  initial  design  can  be  evaluated  using  one  of  the  classification  society 
design  rules.  Appendix  C  and  D  (in  this  case,  API  (Bull  2U)  is  recommended  over  DNV  (RP- 
C202)  due  to  better  predictions  of  failure  pressures).  The  classification  society  design  rules  are  a 
good  check  on  the  analytic  solution  because  of  the  use  of  empirical  data  in  the  derivations  of  the 
equations.  Finally,  the  design  can  be  model  using  MAESTRO™,  which  allows  the  designer 
much  greater  flexibility  in  the  complexity  of  design.  The  MAESTRO™  solution  was  generally 


found  to  be  a  good  first-order  design  approximation  that  predicted  failure  pressures  within  15% 
of  experimental  results.  The  overall  size  and  dimensions  taken  from  MAESTRO™  can  provide 
the  designer  a  good  estimate  of  the  structural  material  required  to  build  the  cylinder  or 
submarine,  which  in  turn  can  be  used  to  approximate  the  structural  weight. 

After  the  initial  design  is  refined  and  judged  to  be  adequate,  a  much  more  rigorous 
analyses  must  be  used  to  ensure  a  safe  design.  These  advanced  analyses  should  include  a  finite 
element  analysis  with  higher  fidelity  (e.g.  ADINA™  or  ABAQUS  ™)  than  die  current  analysis, 
and  other  tools  that  can  provide  a  more  detailed  local  stress  analysis.  These  higher  order  analysis 
tools  can  account  for  material  differences,  geometric  out-of-roundness,  and  actual  construction 
factors  such  as  heat  affected  zones  due  to  welding,  and  bulkhead  effects. 

8.4  Further  Areas  of  Study 

There  are  several  areas  that  require  further  research  to  completely  understand  the  results 
of  the  analysis  conducted  in  this  thesis.  The  area  most  evident  in  need  of  more  study  is  the  large 
difference  between  the  predicted  failure  pressures  determined  using  the  classification  society 
design  rules,  API  (Bull  2U)  and  DNV  (RP-C202),  and  the  actual  failure  pressures.  Due  to  the 
semi-empirical  nature  of  these  design  tools,  it  is  expected  that  larger  test  cylinders  (with  diameter 
to  thickness  ratios  greater  than  300)  would  provide  more  accurate  failure  predictions.  To 
accomplish  this  analysis,  more  test  cylinder  data  would  need  to  be  obtained  with  careful  attention 
paid  to  diameter  to  thickness  ratio,  radius  to  thickness  ratio,  and  minimum  thickness 
specifications. 

There  are  other  classification  societies  that  produce  design  rules  for  stiffened  cylinders 
and  other  geometries.  These  societies  include  NORSOK  (Norway),  Lloyd’s  Register  (United 
Kingdom),  Registro  Italiano  Group  (RINA)  (Italy)  and  several  others.  These  additional  design 


rules  could  be  compared  against  the  existing  test  cylinders.  More  importantly,  the  design  rules 
used  in  [4]  and  the  design  rules  used  in  this  thesis  could  be  compared  using  additional 
experimental  test  cylinder  data. 

Finally,  the  source  code  for  MAESTRO™  was  written  based  on  formulations  from  API 
(Bull  2U,  1987  edition).  The  source  code  should  be  updated  to  reflect  the  current  edition  (2000) 
of  API  (Bull  2U).  MAESTRO™  could  then  be  used  for  comparisons  to  other  finite  element 
programs  (e.g.  ADINA™  or  ABAQUS  ™)  to  determine  the  relative  strengths  and  weaknesses  of 
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Appendix  A:  Analytic  Solution 
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[Analytic  Solution 


NAVSEA  Test  Cylinder  1  .d 


General  Defintions 

ksi  :=  6.894757106Pa  rtog  :=  64.0— 

ft3 


E  :=  3000(ksi 


Young's  Modulus  of  Elasticity 
Inner  radius  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 

Distance  between  bulkheads 
Poison's  ratio  for  Fe/Steel 
Shell  thickness 


kip  :=  lOOOIbf  bload  :=  1 


Lf  =  4.266in 
Lh  =  22.488in 

h  =  0.081in 


Ring  Stiffener  Dimensions 


b :  ^ 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 


flange  thickness  of  ring  stiffener 


r  •=  r.  +  ii.  radius  of  cylinder  to  mid-line  of  shell 
^  2 

D  :=  2-R  diameter  of  cylinder 


=  0.138in 
hw  =  0.57in 


bf  =  Oin 
tf  =  0in 


R  =  8.047in 
D  =  16.095in 


L-Lf-b 


unsupported  shell  length 


L  =  4.128in 


w 


f 


+  tfbf 


c,  :=■ 


***1) 


twhw  +  tfbf 


dist  from  shell  to 
centroid 


Cj  =  0.285in 


o 

1 

* 

A 

II 

£ 

dist  from  centroid  to  end  of  flange 

C2  =  0.285in 

R^g  R  +  .5*h  +  Cj 

radius  to  centroid  of  ring  stiffener 
(external  stiffeners) 

Reg  =  8.373in 

Af:=(twhw+brtf) 

cross-sectional  area  of  ring  stiffener 

Af=0.079in2 

Aeff:=Aff— ^ 

effective  area  of  stiffener 
eqn  [24a]  from  P&S 

Aeff  =  0.076m2 

Aeff 

a  := - 

Lfh 

ratio  of  effective  frame  area  to  shell 
area  eqn  [62]  P&S 

a  =  0.219 

P 


b_ 

k 


ratio  of  faying  width  to  frame  spacing 
eqn  [62]  P&S 


P  =  0.032 


Failure  Modes 

a.  Axisvmmetric  Buckling  (AX) 


Based  on  Poutos  and  Salerno,  1961 


Pcab 


Y  <-  0 

limit  <-  5psi 
test  «->  Opsi 
convert  1  psi 

i «-  o 

while  j  <  10 

1 


n2 

e  <r 

;i  4 


1  +  r 


Aj^r7).-k= 

VrT 

4  cosh  (n  )2  -  cos  (n  2'0)2 

0  cosh  (r\  j  B^sinh  ^  j-e)  cos  (r\  2*0  )  sin  (r\  2-0  ) 


*1  1 

n2 

cosh 

(n  re)  sin  ( 

Tlj-e) 

sinh 

(»1  i-e)  cos  (hi  2-o ) 

n2 

+  - 

’ll 

cosh 

(it  j  G^  sinh 

(n  i0 

J  cos 

(n2e)sin  (»i2  e) 

f4. 


n  i  n2 

cosh  (ti  i*6  )*sin  ^  2*0  )  sinh  ^j  -G^cos  (h2*g) 

J3  _ ^2_ _ _ 

j  _  y2  cosh  (t]  i*G  )*sinh  (t\  |-0  )  cos  2'0  )sin  2-G  ) 

n  i  ii2 

/  \ 


1  V  ^  a 
1 - ,*a 

2j 


A  <_  - V - y  . 

a  +  0  +  (l  -  P  )*F i 
«-  A2|f22  +  F2.F4  (l  -  2  v). 

F2  —  v.F4* 


K1 

k2 


w 


0.91 


v(i) 


i-y2  J 


\f42.(i-v  +  v2).J 


J 


[3 

JT+K1-K2 

break  if  |pc  -  test  |  <  limit 


2-E 
test  <-  pc 
j  <-  j  +  1 


out 


convert 


out  x  <-  j 


out 


f  874.261  ^ 
peab  ={  2  ) 


PcAX  ^Pcab^PS' 
kcAX  =  874.26psi 


b.  Asymmetric  Collapse  (Lobar  Buckling)  (LB) 

Based  on  von  Mises  theory,  with  approximations  by  Windenburg,  1933 

Windenburg  approximation  assumes  number  of  lobes  (n)  is  n*2P  /L 


PcLB := 


Lobar  buckling 
[PNAeqn  19] 


=  623.477psi 


c.  General  Instability  fGh 
Based  on  Kendrick,  1954 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


0 


,Uv). 

Le  :=  \.56yfiTh 


L 

V^h 

'  cosh  (o)  -  cos  m 
^  sinh(0)  +  sin(o)  J 


effective  shell  length,  [P&S  eqn  92] 
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ngi_range  : 


ngi_int  2  ^ 
ngi_int  “  * 
ngi_int 
ngi_int  +  * 
^ngi_int  +  2J 


pcG^ngi_range0j  ^ 

pcG^ngi_range  j  j 

f  148926.73^ 

PcGIrange 

pcG^ngi_range2^ 

PcG^ngi_range3j 
/  \ 

PcGI_range 

5059.37 

1240.24 

1434.37 

v  2157.73  ) 

^pcG^ngi_range4J  j 

1 

PcGI  ■  m^n(PcGI_range) 

mOJIApsi 


m> 


Summary 

Axisymmetric  Buckling  §^^  =1874  2^)51 


Lobar  Buckling 


General  Instability 


PCGI=  1240.24psi 
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Analytic  Solution 
NAVSEA  Test  Cylinder  l.f 


General  Defintions 

ksi  :=  6.894757 106Pa  rtog  :=  64.0— 

ft3 


E  :=  3000(ksi 


Young's  Modulus  of  Elasticity 
Inner  radius  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 

Distance  between  bulkheads 
Poison's  ratio  for  Fe/Steel 
Shell  thickness 


kip:= 


Rina  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


r  =  r.  +  radius  of  cylinder  to  mid-line  of  shell 

*  2 

D  :=  2-R  diameter  of  cylinder 


looabf  bload  := 

Lf  =  2.666in 
1^=42.129111 

h  =  0.263in 

=  0.198in 
hw  =  0.762in 

bf  =  0.763in 
tf  =  0.263in 

R  =  17.328in 
D  =  34.657in 


L  :=  Lf  -  b 


unsupported  shell  length 


L  =  2.468in 


Vhw  .  .  f, 

— +  fbr['-,+7; 

Vhw+,fbf 


dist  from  shell  to 
centroid 


C1 


c2:=hw_cl 

dist  from  centroid  to  end  of  flange 

Rj-g  :=  R  -  .5  h  -  cj 

radius  to  centroid  of  ring  stiffener 
(external  stiffeners) 

c2 


Af  :=  (tw-hw  +  bf  tf)  cross-sectional  area  of  ring  stiffener 


Af 


effective  area  of  stiffener 
eqn  [24a]  from  P&S 


Aeff 

Lfh 


ratio  of  effective  frame  area  to  shell 
area  eqn  [62]  P&S 


P 


ratio  of  faying  width  to  frame  spacing 
eqn  [62]  P&S 


=  0.674in 

=  0.088in 

=  16.523in 

=  0.352in2 

=  0.369in2 

0.526 

0.074 


Failure  Modes 

a.  Axisvmmetric  Buckling  (AX) 


Based  on  Poulos  and  Salerno,  1961 


P  cab 


y  <-  0 
limit  4-  5psi 
test  <-  Opsi 
convert 

i  *-  o 

while  j  £  10 


lpsi 


♦ 

n  2 * 

e  «- 

Ifj  *■ 


-I./TT7 

-i.VTT7 


VR^h 


cosh  | 

- 1 - ; - : - +  - - - 

^  1 

[tj  j-0  )*sin  (t\  2*^  )  sinh  i  *  ©  ) 

^2 
•cos  | 

S  2  e) 

n2  ’i  l 

cosh  | 

^j-8^sinh  cos  [x\  2*8 

)  sin  | 

(n2e) 

M 

cosh  (tj  j- 


^2 


3 

*1  2 

n\ 

.  2  cosh  { 

1  -  v  1 

^  j-8  J* sinh  | 

(n  i'0) 

|  cos  | 

^2-e)-sin  | 

(l2-0) 

f  I  -  — -1-a 

A  _  l  ^ _ 

a  +  p  +  (l  -  P  )’F  j 

f2M'  -2v)(J 

✓  -  x 


0.91 


1-v2  j 


^  +  F42-(i  -  v  +  V2) 
I 


K  2 


Pc 


0.91  _  ^ 
2 ) 


a  fiO 
y  'l  R  J 


Jt+k>"K2 

break  if  |pc  -  test  |  <  limit 


2-E 

test  <-  pc 

j  <-  j  +  1 

out 


convert 


out  j  <-  j 


out 


„  _  |  2.176x  lO3^ 

Pcab  ~  , 

l  2  ) 

PcAX  :=  Pcab0’'Psl 

|cAX  =  21'76.44psi 


b.  Asymmetric  Collapse  (Lobar  Buckling)  {LB} 

Based  on  von  Mises  theory,  with  approximations  by  Windenburg,  1933 

Windenburg  approximation  assumes  number  of  lobes  (n)  is  /L 


2.42*  E- 


PcLB: 


fiO 


3 

r 

G-vf. 

Lobar  buckling 
[PNAeqn  19] 


jcLB=122ll.34psi 


c.  General  instability  (GH 
Based  on  Kendrick,  1954 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 

VrIi 

L.  :=  \.S6jfTh(  ■cosh(e)~cos(e)"|  effective  shell  length,  [P&S  eqn  92] 
V  sinh(0)  +  sin(0)  ) 


Le  =  2.41  in 


m 


m 


16 


d:=h„,+  -  + 


w 


*1 

2 


Ap  Le-h  Aw .  tw-hw  Afj .  tfbj  A j .  Ap  +  Aw  +  A^j 


~J  +  AtfAP 


Ig-AT-d2^ 

Ie  =  0.182in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


» 


* 


General  Instability  Calculation 


Given 

n  >  1 


iigj  :=  Minimiz^pcGI)n) 


=  2.984 


ngi_mt:=round(ngi) 


must  be  integer  value 


General  Instability 
calculation 

[SNAME  eqn  10] 
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girange  * 


ngi_int  * 


ngi_int  +  1  I 
^ngi__int +  2 ) 


I_range  *“ 


PcG^ngi_range0^ 
PcG(n  girange  ^ 
PcG(n  gi_range2^ 
pcG^ngi  range 
PcG^ngi_range4^  j 


^213331.15^ 

11481.17 

PcGIrange  =  441753  Psi 

6173.06 
^  9537.92  J 


pcGi:=min(pcGi  range) 


cGI=  4417  53Psi 


Summary 


Axisymmetric  Buckling 


PcAX 


Lobar  Buckling 


General  Instability 


P:c<JI=4417-53psi 
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[Analytic  Solution 


NAVSEA  Test  Cylinder  2.a 


General  Defintions 


ksi:=  6.894757 10  Pa 


rtog  :=  64.0 


lbf 


CTy  :=  65.5ksi 


v  :=  .3 


ft 


E  :=  3000(ksi  Young's  Modulus  of  Elasticity 
Inner  radius  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 

Distance  between  bulkheads 
Poison's  ratio  for  Fe/Steel 
Shell  thickness 


kip  :=  lOOClbf  bload  := 


Lf=  1.366in 
Lb  =  8.636in 

h  =  0.086in 


Ring  Stiffener  Dimensions 


ggm 

thickness  of  web  of  ring  stiffener 

tyj  -  0.044in 

IBP 

height  of  web  of  ring  stiffener 

hw=0.454in 

b:=tw 

faying  width  of  stiffener  (from  P&S 
for  1  beam  stiffener) 

mm 

breadth  of  flange  of  ring  stiffener 

bf  =  0.399in 

flange  thickness  of  ring  stiffener 

tf  =  0.078in 

R  :=  Rj  +  — 

*  2 

radius  of  cylinder  to  mid-line  of  shell 

R  =  8.41 8in 

D  :=  2R 

diameter  of  cylinder 

D  =  16.836in 

L  :=  Lf  -  b 


unsupported  shell  length 


L=  1.322in 


Vhw  .  (  *f\ 

2  +  tfbf|hw  + 

*w,hw+tfbf 


dist  from  shell  to 
centroid 


Cj  =  0.389in 


c2  :=hw~  C1  dist  from  centroid  to  end  of  flange 

Rp  :=  R  +  .5-h  +  cj  radius  to  centroid  of  ring  stiffener 

(external  stiffeners) 


c2  =  0.065in 


Rpg  =  8.85in 


Af :=  (tw*bw  +  bf  tf)  cross-sectional  area  of  ring  stiffener  Af= 0.05  tin2 


effective  area  of  stiffener 
eqn  [24a]  from  P&S 


veff 


=  0.049in 


Aeff 

o  := - 

Lj-h 


ratio  of  effective  frame  area  to  shell 
area  eqn  [62]  P&S 


a  =  0.415 


ratio  of  faying  width  to  frame  spacing  B  =  0  032 
eqn  [62]  P&S 


Failure  Modes 

a.  Axisvmmetric  Buckling  (AX) 


Based  on  Poulos  and  Salerno,  1961 


Pcab 


-  5psi 
Opsi 


y  <- 

limit 
test 

convert  <-  1  psi 

j  0 

while  j  £  10 


*2 

e  <- 

7i 4 


2  ^  * 


-  fVTTF 

>/rT 

4  cosh  fri'6)2  -  cos  Ql2o)2 _ 

e  cosh  (tj  j-8  )*sinh  ^  1*®)  cos  ^  2*0  )*sin  (tj  2*0 ) 

ni  n2 

cosh  (i)  j-e)-sm  (y\  2*®)  sinh  (ri  j-0  )*cos  (ti  2*0  ) 

n2 _ _ 

cosh  ^T|j*0^  sinh  1  * ®  )  cos  (*1  2*0  )*sin  (t|  2*0 ) 

n  1  n  2 

cosh  p0^*sin  ^2*0)  sinh  ^j-B^-cos  (112*®) 

J  3  _ *2 _ M _ 

J  l  -  v2  00511  (’l  l*e  )  sinh  (11  j-0  )  cos  ^  2*0  )  sin  2’0) 
*11  *1 2 


1  -  •— -l*a 

1± 


\  < _ v _ fL' - 

a  +  0  +  (l  -  p)-F| 

-an— r 


-2S-'1 .  f„’(i 

i-v’; 


A 

V  +  V 


(O 

Mil 


-.--r  1-^) 


1  -  V 


c  fi 

break  if  |pc  -  test  |  <  limit 


2-E 
test  <-  pc 

j  <-  j  +  1 

Pc 

out  < - 

u  convert 
out  4-  j 


# 


(  884.668^ 

Pcab=l  2  ) 

Pc  AX  :=Pcab0-1Psi 
$cAX  =  884.67 psi 


b.  Asymmetric  Collapse  (Lobar  Buckling^  (LB1 

Based  on  von  Mises  theory,  with  approximations  by  Windenburg,  1933 

Windenburg  approximation  assumes  number  of  lobes  (n)  is  W*2P /L 


2.42E- 


PcLB:= 


fiO 

U ) 


ita-ts 


Lobar  buckling 
[PNAeqn  19] 


% 


3113.77psi 


c.  General  instability  (Gil 
Based  on  Kendrick,  1954 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 
4r 


VlTh 

Le  :=  1.56>/rT/  cosh  (6 )- cos  (e)^  effective  shell  length.  [P&S  eqn  92] 
v  sinh(0)  +  sin(e)  ) 


Le  =  1.221  in 
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* 


# 


d:=hw  + 


tf 

2 


Ap:=Le-h  Aw:=twhw  Afl:=tfbf  AT  :=  Ap  +  Aw  +  Afl 


Aw^  . 

— )  +  AffAp 


le  :=  AT  d2-Cj 
Ie  =  0.008in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


General  Instability  Calculation 


n  :=  1 


PcGl(n)  := 


Eh 

R 


2  X2  ^  (  2  l\2 

n  + - l  An  +  X  ) 

LV  2  J 


(n2  -  l)-E  I 
R3Lf 


Given 


n  >  1 

ngi  :=Minimi2^pcGI,n) 


-  3.975  must  be  integer  value 


ngi_int :=  round(ngi) 


General  Instability 
calculation 

[SNAMEeqn  10] 
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gi  range  •" 


PcGI_range :_ 


^ngi_int  “  ^ 
ngi_int  “  1 
ngi_int 
ngi_int  +  1 
i^ngi_int  +  2) 


j  PcG^ngi_range  |  ^ 
jPcG^ngi  range  | 
|PcG^ngi_range2j| 
|PcG^ngi_range3j| 
jPcGl^ngi_range4^|  j 


PcGI 


range ' 


f 20396.34^ 
8552.29 
6391.31 
7626.22 
10294.03 ^ 


psi 


PcGI-  ^KPcGI  range) 


Summary 

Axisymmetric  Buckling 


|^^.67psi 


Lobar  Buckling 


General  Instability 


cGI=  6391.3lpsi 
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Analytic  Solution 
NAVSEA  Test  Cylinder  2.c 


General  Defintions 

ksi  :=  6.894757 106Pa  rtog  :=  64.0— 

ft3 


E  :=  3000(ksi 


Young's  Modulus  of  Elasticity 
Inner  radius  of  cylinder 

Yield  strength 

Length  of  supported  cylinder 

Distance  between  bulkheads 
Poison's  ratio  for  Fe/Steel 
Shell  thickness 


kip  :=  lOOQbf  bload  :=  1 


Lf=3.256in 
1 1  =  115.532m 

h  =  0.337in 


Ring  Stiffener  Dimensions 


b :  ^ 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 


flange  thickness  of  ring  stiffener 


D  :=  2-R 


radius  of  cylinder  to  mid-line  of  shell 
diameter  of  cylinder 


tw  =  0.127in 
hw  =  2.01in 


bf  =  1.552in 
tf  =  0.305in 

R  =  18.883in 
D  =  37.765in 


L  :=  Lf  -  b 


unsupported  shell  length 


L=3.129in 


Ci  := 


*w’^w  ,  (, 

2 

tw‘bw  +  *f kf 


dist  from  shell  to 
centroid 


cj  =  1.757in 


o 

1 

£ 

A 

II 

cS" 

dist  from  centroid  to  end  of  flange 

C2  =  0.253in 

Rcg:=R-.5-h-ci 

radius  to  centroid  of  ring  stiffener 
(internal  stiffeners) 

Rpg  =  16.957in 

Af:=  (twhw+  bftf) 

cross-sectional  area  of  ring  stiffener 

Af=0.729in2 

Aeff:=Aff— ^ 

effective  area  of  stiffener 
eqn  [24a]  from  P&S 

Aeff=0-811in2 

Aeff 

a  := - 

Lfh 

ratio  of  effective  frame  area  to  shell 
area  eqn  [62]  P&S 

a  =  0.739 

P:=w 

% 

ratio  of  faying  width  to  frame  spacing 
eqn  [62]  P&S 

P  =  0.039 

Failure  Modes 

a.  Axisvmmetric  Buckling  (AX) 


Based  on  Poulos  and  Salerno,  1961 


Pcab  := 


y  <-  0 
limit  + 
test  <- 
convert 

j  <-  o 

while 


^1 

^2 


-  5psi 
Opsi 
<-  lpsi 

j  <M0 

-  r 


F4 


4  cosh  (ti  j 

•0  )2  -  cos  (tj  2*0  )2 

0  cosh  j -  0  ^  sinh  (r)  j-fl)  cos  (r\  2*0  )• 

sin  ^2*0) 

cosh  ^  j-0  ^-sin  (ti  2*e 

*1 

|  sinh  (tj  j*0  y«>s 

2 

(1l2e) 

*12 

*1  1 

cosh  ^T|  j  * 0  ^ ‘ sinh  j-( 

)  ^  cos  (“H  2 *  ®  )’s*n 

(n2e) 

*1  1 

cosh  ^TJ  j-0  J- 

*12 

sin(r|2*0)  sinh  ^tij-O^  cos 

13  2 

^1 

-Jt  -  V2  00811  (’ll  6) 

sinh  (n  r  0  )  cos  ^ 

r]  2'0  )sin  (n  2*0  ) 

^1  ^2 

f  1  -  —  Yet 

A  <_  _ i  _ 

a  +  P  +  0  -  P)  F, 

^,0- 

(3}  .  (c  K  I  0.91  ^ 

*2^  It  iaI  F2  -  vF4-  I - 7 


PC  +• 


«m . -  js; 

•.ft) 

il+K«-K2 

test  |  ^  limit 


break  if  |pc  f 


y  ■> 

2-E 

test  <-  pc 
j  +-  j  +  1 


*  V 

out  < - 

u  convert 

out  .  j 


out 


„  _  I  4.137  x  103  ^ 

Pcab  “  . 

V  2  J 

PcAX  :=Pcab0-,Psi 
|cAX  =  4137.25  psi 


b.  Asymmetric  Collapse  (Lobar  Buckling)  (LB) 

Based  on  von  Mises  theory,  with  approximations  by  Windenburg,  1933 

Windenburg  approximation  assumes  number  of  lobes  (n)  is  n*2 P  /L 


2.42E- 


PcLB  := 


3 

r 

.GM*)’ 

Lobar  buckling 
[PNAeqn  19] 


||cLB  =14M26psi 


c.  General  Instability  (Gl) 
Based  on  Kendrick,  1954 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


Le:=  1.56 ypThl cosh  (e)  ~  cos  (e)  ^ 
V  sinh(o)  +  sin(e)  ) 


effective  shell  length,  [P&S  eqn  92] 


Le  =  3.03  in 
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gi_range  *“ 


^ngi_int  ^ 
ngi_int  “  1 
ngi_int 
ngi_int +  1 

KiJnt +  2  J 


PcGI 


l_range  *" 


PcGl^ngi_range0^|  ^ 
PcG^ngi_range  j^| 
PcG^ngi_range2j| 
PcG^ngi_range3j| 
PcGl^ngi  range  ^  |  j 


PcGI_  range 


( 619329.53^ 
176802.28 
8506.69 
20994.71 
39274.39 ) 


psi 


PcGI :  m^n(PcGI_range) 


icGI=  8506.69jpsi 


Summary 

Axisymmetric  Buckling 


Lobar  Buckling 


14528.2^si 


General  Instability 


|e>cGi  ~  8506.69psi 


90 


Appendix  B:  Numerical  Solution 


Page  Intentionally  Left  Blank 


Numerical  Solution 
(MAESTRO,  version  8.5) 

NAVSEA  Test  Cylinder  l.d 

General  Defintions 


E  :=  30000 ksi 


v  :=  .3 


Rj  :=  8.00'3n 


D0:_  2(Ri+t) 


ksi  :=  6.89475710  Pa 


Young’s  Modulus  of  Elasticity 
Yield  Strength 

Poison’s  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 

diameter  to  outside  of  shell 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


D0  =  16.176in 


radius  to  centerline  of  shell 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

=  0.138in 

height  of  web  of  ring  stiffener 

hw  =  0.57in 

j* 

ii 

faying  width  of  stiffener  (from  P&S 
for  1  beam  stiffener) 

~  ; 

breadth  of  flange  of  ring  stiffener 

bf  =  Oin 

flange  thickness  of  ring  stiffener 

,3 

o 

II 

u-* 

MAESTRO  Modeler  representation  of  cylinder  l.d 
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HYDROSTATIC  TEST  -  CYL  l.D 


VERSION  OF  DATA  SET  IS  8.1.1 

JOB  TYPE:  STRUCTURAL  ANALYSIS 

DEFLECTIONS  TO  BE  SAVED  ON  FILE  CYL  1D-2.DEF 

STRUCTURE  PARAMETERS: 

********************* 


PLOT  LEVEL  =  2 

TRANSVERSE  SYMMETRY  INDICATOR  =  2 

LEVEL  OF  OUTPUT  REGARDING  F.  E.  MODEL  =  1 

DEFAULT  EVALUATION  LEVEL  =  3 

STATION  SPACING  FOR  PRINTING  SUMMED 

VERTICAL  LOADS  =  1st  sectn. 

FIRST  SUBSTRUCTURE  IN  NODE  RENUMBERING:  1 

FIRST  MODULE  IN  NODE  RENUMBERING:  1 

GLOBAL  X  VALUE  FOR  STATION  0:  Lowest  value 

IF  1,  SUPPRESS  OUT-OF-PLANE  DEFLECTION 

OF  UNSTIFFENED  COMPOUND  NODES:  0 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR 
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Numerical  Solution 
(MAESTRO,  version  8.5) 

NAVSEA  Test  Cylinder  l.f 


General  Defintions 


E  :=  30000 ksi 


17.191b 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison's  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 


ksi:=6.894757106Pa 


diameter  to  outside  of  shell  D0  =  34.92in 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


radius  to  centerline  of  shell 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw  =  0.198in 
hw  =  0.762in 


bf  =  0.763in 
tf  =0.263in 


MAESTRO  Modeler  representation  of  cylinder  l.f 
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HYDROSTATIC  TEST  -  CYL  l.F 

VERSION  OF  DATA  SET  IS  8.1.1 

JOB  TYPE:  STRUCTURAL  ANALYSIS 

DEFLECTIONS  TO  BE  SAVED  ON  FILE  CYL  1F.DEF 

STRUCTURE  PARAMETERS :  # 

********************* 


PLOT  LEVEL  = 

TRANSVERSE  SYMMETRY  INDICATOR  = 

LEVEL  OF  OUTPUT  REGARDING  F.  E.  MODEL  as 
DEFAULT  EVALUATION  LEVEL  = 

STATION  SPACING  FOR  PRINTING  SUMMED 
VERTICAL  LOADS  = 

FIRST  SUBSTRUCTURE  IN  NODE  RENUMBERING: 
FIRST  MODULE  IN  NODE  RENUMBERING: 

GLOBAL  X  VALUE  FOR  STATION  0: 

IF  1,  SUPPRESS  OUT - OF - PLANE  DEFLECTION 
OF  UNSTIFFENED  COMPOUND  NODES: 
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REFERENCE  COORDINATES  OF  STRUCTURE  ORIGIN  AND  OCEAN  SURFACE 
KEYWORD  XREF .  YREF.  ZREF .  OCEAN  SURFACE  w 
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Numerical  Solution 
(MAESTRO,  version  8.5) 


NAVSEA  Test  Cylinder  2.a 


General  Defintions 

E  :=  30000 ksi 


v  :=  .3 


Young's  Modulus  of  Elasticity 
Yield  Strength 


ksi  :=  6.89475710  Pa 


Poison's  ratio  for  Fe/Steel 
thickness  of  shell 


Dg  =  16.922in 


r.  :=  8.375n  Inner  radius  of  cylinder 

D0  :=  2-(Rj  +  t)  diameter  to  outside  of  shell 

ring  spacing  (frame  center  to  frame  center) 
length  of  cylinder  between  bulkheads  or  lines  of  support 


radius  to  centerline  of  shell 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw  =  0.044in 
hw=0.454in 


bf=0.399in 
tf  =  0.078in 
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MAESTRO  Modeler  representation  of  cylinder  2.a 


MAESTRO  Version  8.5.  0  12-APR— 2003 

HYDROSTATIC  TEST  -  CYL  2 .A 

VERSION  OF  DATA  SET  IS  8 . 1 . 1 

JOB  TYPE:  STRUCTURAL  ANALYSIS 

DEFLECTIONS  TO  BE  SAVED  ON  FILE  CYL  2A.DEF 


PAGE 


STRUCTURE  PARAMETERS: 
********************* 


1 


PLOT  LEVEL  = 

TRANSVERSE  SYMMETRY  INDICATOR  = 

LEVEL  OF  OUTPUT  REGARDING  F.  E.  MODEL  = 
DEFAULT  EVALUATION  LEVEL  = 

STATION  SPACING  FOR  PRINTING  SUMMED 
VERTICAL  LOADS  = 

FIRST  SUBSTRUCTURE  IN  NODE  RENUMBERING: 
FIRST  MODULE  IN  NODE  RENUMBERING: 

GLOBAL  X  VALUE  FOR  STATION  0: 

IF  1,  SUPPRESS  OUT-OF-PLANE  DEFLECTION 
OF  UNSTIFFENED  COMPOUND  NODES: 


2 

2 

1 

3 


1st  sectn. 

1 

1 

Lowest  value 


REFERENCE  COORDINATES  OF  STRUCTURE  ORIGIN  AND  OCEAN  SURFACE 
KEYWORD  XREF.  YREF.  ZREF .  OCEAN  SURFACE 

REFERENCE  0.000  0.000  0.000  0.000 
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DEFINITION  OF  STRAKES 
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LOADCASE:  HYDROSTATIC  PRESSURE  -  797  PSI 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  3  OF  SUBSTR 
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Numerical  Solution 
(MAESTRO,  version  8.5) 

NAVSEA  Test  Cylinder  2.c 


General  Defintions 

E  :=  30000 ksi  Young's  Modulus  of  Elasticity  :=  6-89475710  Pa 

Yield  Strength 


v  ;=  3  Poison's  ratio  for  Fe/Steel 

HUjg^ll  thickness  of  shell 

Rj  :=  18.714n  Inner  radius  of  cylinder 

D0  :=  2-^  +  t)  diameter  to  outside  of  shell  D0  =  38.102m 

ring  spacing  (frame  center  to  frame  center) 
length  of  cylinder  between  bulkheads  or  lines  of  support 


radius  to  centerline  of  shell 


R:= 


_2  _1 
2  2 


Ring  Stiffener  Dimensions 


b:=t 


w 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw=  0.127in 


hw  =  2.01in 


bf  =  1.552in 
tf  =  0.305in 
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MAESTRO  Modeler  representation  of  cylinder  2.c 


Y 


X 


MAESTRO  Version  8.5.  0  12-APR-2003  PAGE  1 

HYDROSTATIC  PRESSURE  -  CYL  2C 


VERSION  OF  DATA  SET  IS  8.1.1 

JOB  TYPE:  STRUCTURAL  ANALYSIS 

DEFLECTIONS  TO  BE  SAVED  ON  FILE  CYL  2C.DEF 

STRUCTURE  PARAMETERS: 

********************* 


PLOT  LEVEL  = 

TRANSVERSE  SYMMETRY  INDICATOR  = 

LEVEL  OF  OUTPUT  REGARDING  F.  E.  MODEL  = 
DEFAULT  EVALUATION  LEVEL  = 

STATION  SPACING  FOR  PRINTING  SUMMED 
VERTICAL  LOADS  = 

FIRST  SUBSTRUCTURE  IN  NODE  RENUMBERING: 
FIRST  MODULE  IN  NODE  RENUMBERING: 

GLOBAL  X  VALUE  FOR  STATION  0: 

IF  1,  SUPPRESS  OUT-OF-PLANE  DEFLECTION 
OF  UNSTIFFENED  COMPOUND  NODES: 


2 

2 

1 

3 


1st  sectn. 

1 

1 

Lowest  value 


REFERENCE  COORDINATES  OF  STRUCTURE  ORIGIN  AND  OCEAN  SURFACE 
KEYWORD  XREF.  YREF.  ZREF.  OCEAN  SURFACE 

REFERENCE  0.000  0.000  0.000  0.000 
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IV (A)  -  DEFINITION  OF  STRAKES 
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Appendix  C:  API  (Bull  2U)  Solution 
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Classification  Society  Solution 
(API  Bulletin  2U) 


NAVSEA  Test  Cylinder  l.d 


General  Pefintions 


ksi:=  6.89475710  Pa 


rtog  :=  64.0—  kip  :=  lOOdbf  bload  :=  1  — 


E  :=  30000 ksi 


v  :=  .3 


r 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison's  ratio  for  Fe/Steel 
thickness  of  shell 

gllllllll  Inner  radius  of  cylinder 
Dq  :=  2-(Rj  +  t)  diameter  to  outside  of  shell 


kpa  := 


in 

kg 


2 

sec  m 


D0=  16.176in 


ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw  =  0.138in 
hw  =  0.57in 


bf  =  0  in 
tf  =  Oin 
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*°-  2 
R:=Ro-i 


Vh 


w 


radius  to  outside  of  shell 
radius  to  centerline  of  shell 


c,  :=■ 


(  •ft 

VNv  +  V  bf 


dist  from  shell  to 
centroid 


Cj  =  0.285in 


c2:=hw-cl 


dist  from  centroid  to  end  of  flange 


distance  from  centerline  of  shell  to  centroid 
of  ring  stiffener  (bositive  outward ) 


Rj.  :=  R  +  Zf  radius  to  centroid  of  ring  stiffener 
Aj  :=  +  bftf)  cross-sectional  area  of  ring  stiffener 


C2  =  0.285in 
Zj.  =  0.325in 


A,  =  0.079in2 


L^Lj-b 


unsupported  shell  length 


L  =  4.128in 


1  3  (  hw^2  j  , 

*r:=  +  W^cl  ~~J  +  72  bf tf  +  bf * 


■  tfV 
°2+l) 


I,.  =  2. 13  X  10  3  in4  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


v  sinh(o)  +  sin(e)  J 


effective  shell  length,  [P&S  eqn  92] 
Le=  1.255in 
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Afl:=tfbf 


AT :=  \  +  Aw  +  Afl 


* 


Aw Sv'b 


w 


O 

~Tj  +  AflAp 


yf 


2 

2 


tf  +  d-C^ 


Ie  •- AT-d2-Cj 
Ie  =  0.007in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


«D 


Rc:=R_  2 +yP 


radius  to  centroidal  axis  of  combined  ring 
stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


Rj,  =  8.181in 


* 
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'Failure  Modes 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 
4.1 .1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

h_ 

D:=2R 


Mx:=- 


Mx=  5.284 


D  =  16.095in 


c  := 


2.64  if  Mx£  1.5 


3.13 


W 


0.42 


if  1.5<  Mv<  15 


1.0  if  MXS>  15 


axL:_ 


0.207  if  —  ^  1242 
t 


(169c) 


195  + 


Kt)  ' 


if  —  <  1242 


c  =  1.556 


«xL=  0893 


Cx:=  0.605  for  D/t  >  300 


FxeLr-  axL-(-x'2'E  D 


FxeLr=  U25x  «TPa 
b  -  Inelastic  Buckling  Stresses 


rxcLra' 


233F„ 


166+0.5^yJ 


if  —  <  600 
t 


(0.5Fy)  if  y  k  600 


pxcLr 


:=  min((  F, 


xeLr  rxcLra, 

.8, 


)) 


FxcLr=4843x  10  Pa 
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Choose: 


k  =  0.5  k  =  0  for  radial  pressure  and  0.5  for  hydrostatic  pressure 


a  -  Elastic  Buckling  Stresses 

Am:=Mx-1.17+  1.068c 

(“J 


Am=  4.648 


Cp  := 


Cp Cp 


cp  =  0.047 


#> 


PeL:= 


1.18  .  VDj 


Am  +  0-5 


3-*f±y 


if  Mx>  1.5  a  Am<  2.5 


if  2.5<Am<0.104^ 


|^6.688Cp  1061-e  ^~j  if  0.208  <Cp<  2.85 
2.2E j  if  Cp  £  2.85 


a6L:=  8 


A  :=  A, 


^Rf 

a; 


imperfection  factor,  normally  0.8  for  fabricated  cyl 


Le  :=  1.567^  +  V 


A  =  0.073in2 


Le  =  1.397in 


section  11.3 


V  := 


1.0  if  Mx<  1.26 


(l.58-0.4flvl^  if  1.26<MX<3.42 


0  if  Mx>3.42 


1  -  0.3k 


e  :=  * 


1  + 


v* 
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KflL:=  1.0  if  MXS  3.42 
1  -  e\|/ 


a0L‘PeL'®O 


FreLr=  48057ksi 


c  -  Inelastic  Buckling  Stresses 


Ac  =  0.601 


r\  :=  1 1  if  Ac  £  0.55 


0.45  ^ 

- +  0.18  if  0.55  <  Ac  £1.6 

Ac  ) 


1  +  1.15Ac 


if  1.6  <  Ac  <6.25 


—  if  Ac  £  6.25 
Ac 


FrcLr:  ^  FrcLr 


ti  =  0.929 


FrcLr=  4-465x  104  psi 


d  -  Failure  Pressures 

PcLr:=TVa9L'PeL 


•iPcLr=447-)6^)si 
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4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 
4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


Arb  :== ' 


axL:= 


0.207  if  —  >  1242 
t 

(169)  D 

- — L -  if  —  <  1242 


195  + 


axG:= 


0.72  if  Arb  >  0.2 


[(3 -6  ~  5  &axl)' Arb  +  axi]  if  0  06 <  Arl 


axL  if  Arb  -  0  06 


FxeGr:=axG°-605E^j(1  +  Arb) 


FxeGr=  1  005*  10  Pa 


c  -  Inelastic  Buckling  Stresses 


Ac : 


rxeGr 


r\  := 


1  if  Ac  <  0.55 

—  +  0.18)  if  0.55 <  Ac  <  1.6 
Ac  J 


( 


1.31 


1+  l.lSAc 

—  if  Ac  >  6.25 
Ac 


if  1.6<  Ac  <6.25 


^xcGr  ^'^xeGr 
FxcGr=  4-253  x  108Pa 


<0.2 
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a  -  Buckling  Stresses  With  or  Without  End  Pressure 


v= 


(l.l if  Mx>  1.56 
l*  if  Mx^  1.56 


-  T  A  72  Le‘  Le  * 

U  •=  I*.  +  AZ - + - 

er  '  r  r  Af+Le-t  12 


-  JtR 
^:=7" 

h 


PeG<n):=- 


E(iK 


EIcr  (»2  -  ') 


O^kV-.X.2^)2  w** 


n  :=6 


Given 
2£n  <  15 


nj  :=Minimiz^peQ,n) 
nj  =  3.241 


i 


K0G:=(l  ~  0.3k)- 


V 


Ar+Let 


a0G:=  8 


imperfection  factor,  normally  0.8  for  fabricated  cyl 


FreGr:=a0G 


Pednl)^0 


*0G 


FreGr=3-37*  10  Pa 


m 


# 
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c  -  Inelastic  Buckling  Stresses 

FreGr 
Ac - 

Fy 

1  if  Ac  <0.55 

(—  +  0.18!  if  0.55 <  Ac  <;  1.6 

l  Ac  ) 

1.31 

-  if  1 .6  <  Ac  <  6.25 

1  +  l.lSAc 

—  if  Ac  >  6.25 
Ac 

FrcGr=  ^‘^reGr 

FrcGr=309x  1()8pa 
d  -  Failure  Pressures 

PcGr:=TVa0G-PeG(nl) 


SUMMARY 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 

4.1.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

FxeLr=  1125x  109  Pa 
b  -  Inelastic  Buckling  Stresses 

FxcLr=  4843x  108  Pa 

4.1.2  -  External  Pressure 
a  -  Elastic  Buckling  Stresses 

FreLr=3313x  1()8pa 
c  -  Inelastic  Buckling  Stresses 

FrcLr=3079x  1()8pa 
d  -  Failure  Pressures 

pcLr=3083x  106  Pa 

4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 

4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

FxeGr=  1005x  1q9  Pa  FxeGr=  1  005x  1()9pa 

c  -  Inelastic  Buckling  Stresses 

FxcGr=  4253  x  108  pa  FxcGr=  4253  x  10®  Pa 

4.2.2  -  External  Pressure 

♦  a  -  Buckling  Stresses  With  or  Without  End  Pressure 

FreGr=  3368x  1()8pa  FreGr=  3368x  1()8pa 

c  -  Inelastic  Buckling  Stresses 

FrcGr=  3088x  1()8pa  FrcGr=  3  088x  1()8pa 

d  -  Failure  Pressures 

pcGr=^174x  *®6pa  |PcGr=  ^95.44psi 


FxeLr=  1125x  1<)9pa 
FxcLr=4843x  1(,8pa 

FreLr=3313x  1()8pa 
FrcLr=  3079x  108pa 
PcLr  =  447.l6^)si 


m 


» 


m 
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Classification  Society  Solution 
(API  Bulletin  2U) 

NAVSEA  Test  Cylinder  1  .f 


General  Defintions 


6  Ibi 

ksi  :=  6.89475710  Pa  rtog  :=  64.0 —  kip  :=  lOOabf 

ft3 


E  :=  30000 ksi 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison’s  ratio  for  Fe/Steel 
thickness  of  shell 

-  Inner  radius  of  cylinder 


bload 


sec  m 


diameter  to  outside  of  shell  Dq  =  34.92in 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


b:=tw 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 


flange  thickness  of  ring  stiffener 


tw  =  0.198in 
hw  =  0.762in 


bf  =  0.763in 
tf  =  0.263in 


R.-  — 

*0-  2 

R:=Ro-f 


radius  to  outside  of  shell 


radius  to  centerline  of  shell 


Vbw  ( 

2  +  V  f  +  2  J  dist  from  shell  to 


Vhw+tfbf 


centroid 


Cj  =  0.674in 


c2:=hw-ci 


dist  from  centroid  to  end  of  flange 


C2  =  0.088in 


1)  distance  from  centerline  of  shell  to  centroid 

of  ring  stiffener  (positive  outward )  zr =  _0-805in 


Rj.  :=  R  +  Zf  radius  to  centroid  of  ring  stiffener 


Rj.  =  16.523in 


Ar:=(tw  liw  +  bftf)  cross-sectional  area  of  ring  stiffener  Ar  =  0.352in2 


L  :=  Lj.  -  b 


unsupported  shell  length 


L  =  2.468in 


f  _  1  3  (  K')2  1  3  .  (  *fV 

fr  - 12'Vhw  +  Vhw'^cl  "  ~ j  +7^'bftf  +bftf^c2+Jj 


Ij.=  0.03  tin 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


0  :=V3*\1  -  v 


2  L 


Le  :=  1 .5 6ypTt  ( .cosh(9)  cos(e)"\  effective  shell  length,  [P&S  eqn  92] 
V  sinh(0)  +  sin(0)  J 


t  *f 
d  :=  hw  +  —  +  — 
w  2  2 


Le  =  2.41  in 
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Le-t 


Aw :  ^w'^1 


w 


Afl:=tfl 


AT :=  Ap  +  Aw  +  Afl 


m 


Awl 

"Tj  +  AflAp 


Aw 

1t  +  ap/ 

Ax 


°2) 


yf 


-i-tf+dC2 


Ie  :=  AT-d2Cj 
Ie  =  0.182in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


V-K+j-yp 


radius  to  centroidal  axis  of  combined  ring  =  I7.04lin 

stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


# 


# 
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Failure  Modes 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 
4.1.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


D  :=2R 


c  := 


2.64  if  Mx£  1.5 
3.13 

-  if  1.5<  Mv<  15 

0.42  x 


w 


1.0  if  Mx£  15 


axL:= 


0.207  if  —  2>  1242 
t 


(169  c) 


195  + 


Kt)  1 


if  —  <  1242 


Cx:=  0.605  for  D/t  >  300 

FxeLr:=  axL'^x'2F' ^ 

FxeLr=  3248  x  109Pa 

b  -  Inelastic  Buckling  Stresses 


rxcLra' 


233F„ 


166+0.5^J 
(0.5Fy)  if  H 


if  -  <  600 
t 


£  600 


FxcLr  m^(^xeLr  FxcLra)) 
FxcLr=6824x  1<)8pa 


Mx=  1.249 
D  =  34.657in 


c  =  2.64 


axL 


=  1.71 


# 


* 


♦ 
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Choose: 


k  =  0.5  k  =  0  for  radial  pressure  and  0.5  for  hydrostatic  pressure 


m 


* 


# 


a  -  Elastic  Buckling  Stresses 


Am:=Mx-U7+  1068c 


—  =  131.776 
t 


Cp  := 


(?) 


Am=  0.613 


Cp Cp 


Cp  =  0.009 


PeL:= 


_ E.(±) 

1.18  ^  ^  vdJ 


Ajn  +0.5 


a„  Id; 


if  Mx>  1.5  a  Am<2.5 


6.688  CL 


if  2.5  <  Am<  0.104J  --j 


R) 


,0Mi)3 


if  0.208  <  Cp  <  2.85 


2.2-E- 


IV 


if  Cp  >  2.85 


Mx  is  too 
small, 
not  in 
range 


“Value  not  calculated 
(Mx  outside  limits) 


a0L :=  0-8 


Iv 

Le  :=  1.56>/iu  +  t w 


imperfection  factor,  normally  0.8  for  fabricated  cyl 


A  =  0.387in2 


Le  =  3.528in 


section  11.3 


\y:= 


1.0  if  Mx<  1.26 


(1.58-0.46VIJ  if  1.26  <MX<  3.42 


0  if  Mx>  3.42 


1  -  0.3k 


e  :=  ■ 


1  + 
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1.0  if  Mx^  3.42 


1  -  ey 


„  “GL'PeL  ^O  Tr 
hreLr:~  “  ^0L 

FreLr=,Pa 


c  -  Inelastic  Buckling  Stresses 


Ac  := 


rreLr 


r\  := 


1  if  Ac  ^  0.55 
0.45  ^ 

- +  0.18  if  0.55  <  Ac  £1.6 

Ac  ) 

1.31 


(' 


1  +  l.lSAc 

—  if  Ac  £  6.25 
Ac 


if  1.6  <  Ac  <  6.25 


FrcLr-  4  preLr 

FrcLr=  1  Pa 

d  -  Failure  Pressures 

PcLr:=Tla0L'PeL 
P,cLr~ 1  Psi 


"Value  not  calculated 
(Mx  outside  limits) 
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4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 
4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

Ar 


Arb:= 


axL:= 


0.207  if  —  >  1242 
t 

<^—  if  »<!242 


195  + 


K”)  * 


axG:= 


0.72  if  Arb  2>  0.2 

[(3-6  ~  50-axl)-Arb  +  «xl]  if  0  06  <  Ar1 
axL  if  Arb  -  0  06 


FxeGr:=axG°-605E‘ 


FxeGr=1676x  10  Pa 
c  -  Inelastic  Buckling  Stresses 


Ac  := 


rxeGr 


ri  - 


1  if  Ac  <  0.55 

—  +  0.18)  if  0.55 <  Ac  <  1.6 
Ac  ) 

1.31 


( 


1  +  1.15Ac 


if  1.6  <  Ac  <6.25 


—  if  Ac  >  6.25 
Ac 


*xcGr  ^‘^xeGr 
FxcGr=572x  1()8pa 


<1 


<0.2 
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a  -  Buckling  Stresses  With  or  Without  End  Pressure 


v= 


(l.lVlM  +  v)  if  Mx>  1.56 
Lj.  ifMx£1.56 

,  Vt  L„t3 
Ier  •'=  ^  +  Ar-Zr 


Ar+Le-t  12 


.  nR 


PeG^n)  := 


KiK 


BIer(n2-l) 


n  :=  6 

Given 
2  ^  n  <  15 


:=  Minimiz^p^n) 
n2  =  2.963 


V 


Ar+Let 


KeG:=(l-0.3k). 

a0G:=°.8 

peG^nlJ 
*reGr:=  a0G - ~ - ^GG 

FreGr=9-49*  l<>8Pa 
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c  -  Inelastic  Buckling  Stresses 


Ac  :=- 


rreGr 


1  if  Ac  <0.55 

(—  +  0.18)  if  0.55<Ac<  l.f 
l  Ac  ) 


1.31 


1+  l.lSAc 

—  if  Ac  >  6.25 
Ac 


if  1.6  <  Ac  <  6.25 


FrcGr:=1VFreGr 
FrcGr=4-76x  1()8pa 


d  -  Failure  Pressures 

PcGr:=1laeGPeG(nl) 


SUMMARY 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 


4.1.1  -  Axial  Compression  or  Bending 

a  -  Elastic  Buckling  Stresses 

FxeLr=  3248  x  109Pa 

FxeLr=  3248x  1()9pa 

b  •  Inelastic  Buckling  Stresses 

FxcLr=  6824  x  1()8pa 

FxcLr=6-824x  1()8pa 

4.1.2  -  External  Pressure 

a  -  Elastic  Buckling  Stresses 

FreLr=»Pa 

FreLr= ' Pa 

c  -  Inelastic  Buckling  Stresses 

FrcLr=»Pa 

FrcLr=  1  Pa 

d  -  Failure  Pressures 

Pf-T  r  =  *  Pa 

**Vajue  not  calculated 

(Mx  outside  limits) 

General  Instability  of  Ring  Stiffened  Cylinders 

4.2.1  -  Axial  Compression  or  Bending 

a  -  Elastic  Buckling  Stresses 
FxeGr=  1.676x  10?  Pa 

FxeGr=  1676x  1()9pa 

c  -  Inelastic  Buckling  Stresses 

FxcGr=5-72x  lo8pa 

FxcGr=5-72x  1()8pa 

4.2.2  -  External  Pressure 

a  -  Buckling  Stresses  With  or  Without  End  Pressure 

FreGr=9-487*  1()8pa 

FreGr=9-487x  1()8pa 

c  -  Inelastic  Buckling  Stresses 

FrcGr=  4-764x  1()8pa 

FrcGr=4-764x  1()8pa 

d  -  Failure  Pressures 

7 

pcGr=  1.267x  10  Pa 

f*cGr=  1 838.2Cpsi 
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Classification  Society  Solution 
(API  Bulletin  2U) 


NAVSEA  Test  Cylinder  2.a 


General  Defintions 


ksi:=  6.89475710  Pa 


rtog  :=  64.0 


Ibf 


v  :=  .3 


r 


kip  :=  lOOQbf  bload:=l^ 

in 


E  :=  30000ksi  Young's  Modulus  of  Elasticity 

Yield  Strength 


kpa  := 


kg 


2 

sec  -m 


Poison's  ratio  for  Fe/Steel 
thickness  of  shell 


Inner  radius  of  cylinder 
Dq  :=  2-^  +  t)  diameter  to  outside  of  shell 


Dq  =  16.922in 


ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


b:=tw 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


t^,  =  0.044in 
hw  =  0.454in 


bf=0.399in 
tf  =  0.078in 


R°  "  2 

R:=Ro-i 


radius  to  outside  of  shell 


radius  to  centerline  of  shell 


tw-hw 


+  tf-bf-  h,„  +  — 


+  lfDf^nw 

W  +  tfbf 


2  )  dist  from  shell  to 


centroid 


Ci  =  0.389in 


c2:=hw-cl 


dist  from  centroid  to  end  of  flange 


C2  =  0.065in 


Zr  :=  c  j  +  -  distance  from  centerline  of  shell  to  centroid 
2  of  ring  stiffener  (positive  outward ) 

R,.  :=  R  +  Zr  radius  to  centroid  of  ring  stiffener 
Ar  :=  (tw  hw  +  bftf)  cross-sectional  area  of  ring  stiffener 


Zr  =  0.432in 


Af  =  0.05  tin2 


L  :=  bj-  -  b 


unsupported  shell  length 


L=  1.322in 


1  3  (  bw^  1  3  (  tf'] 

*r :=  ^'V^w  +  ~  ~J  +  "j2  bf  tf  +  ^^^2  +  ) 


Ij.  =  l.219x  10  3  in4  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


VR-t 


Le  :=  1.56jRrt  f  — h..(ej  cos(e)^l  effective  shell  length,  [P&S  eqn  92] 
e  ^  sinh(e)  +  sin(e)  ) 


Le  =  1.221in 


t 

d:=hw+-  +  - 
w  2  2 
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Ap  Le-t 


Aw Sv'^w 


Afl:=tfbf 


AT:=Ap  +  Aw+  Afl 


+  AflAp 


yP:=2t  +  d'(1_C2) 


yf 


-~tf+  d-C2 


Ie:=Axd2C1 
Ie  =  0.008in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


radius  to  centroidal  axis  of  combined  ring  =  8.559in 

stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


♦ 
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Failure  Modes 


4.1 


Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 
4.1 .1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

Mx:=-pz 

x  yfTt 

D:=2R 


c  := 


2.64  if  Mx^  1.5 


3.13 


N 


0.42 


if  1.5<  Mv<  15 


a 


xL-: 


1.0  if  MXS  15 


0.207  if  —  £  1242 
t 


( 1 69-  c) 


195  + 


Kt)  1 


if  -  <  1242 


:=  0.605  for  D/t  >  300 
pxeLr:=  q 

FxeLr=  1886x  1()9pa 


b  -  Inelastic  Buckling  Stresses 
233F. 


‘"xcLra" 


y  D 
- -  if -<600 


166+0.5^  * 

(o.5Fy)  if  ^  ^  600 


pxcLr  m*n((pxeLr  pxcLra)) 
FxcLr=  3.984x  108Pa 


Mx=  1.607 
D  =  16.836in 


c  =  2.564 


axL=  1.479 


e 
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Choose: 

k  =  0.5  k  =  0  for  radial  pressure  and  0.5  for  hydrostatic  pressure 


a  -  Elastic  Buckling  Stresses 

Ajn^Mx-  1.17+  1.068k 

(2AJ 


Am=  0.971 


Cp:= 


Cp:=Cp 


) 


Cp" 


=  9.9  x  10 


-3 


PeL:= 


5.08 


ab"8  +  o,  W 


if  Mx  >  1.5  a  Am  <  2.5 


M.s.E.fi'l 


Am  VDj 


if  2.5  <  Am  < 


0104(f) 


“88s‘1061'B(i)  j 


if  0.208  <  Cp  <  2.85 


22H)  j 


if  Cp  £  2.85 


a0L :=  0  8  imperfection  factor,  normally  0.8  for  fabricated  cyl 


A:=Ar- 


a; 


A  =  0.046in2 


section  11.3 


Le 1.56^R-t  + 


Le  =  1.37in 


v|/:= 


1.0  if  Mx<  1.26 


(l.58-0.4<Mx)  if  1.26 <MX< 3.42 


0  if  Mx>3.42 


e  := 


1  -  0.3k 


1  + 


V‘ 
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%/= 


1.0  if  Mx£  3.42 


1  -  £*\|f 


„  “GL'PeL1^  Tr 
FreLr:= - ~ - 

FreLr=  1  172x  109Pa 


c  -  Inelastic  Buckling  Stresses 


Ac  := 


‘reLr 


r\  := 


1  if  Ac  <0.55 
( 0.45  ^ 

- +  0.18  if  0.55  <  Ac  £  1.6 

l  Ac  ) 


1.31 


1  +  l.lSAc 

—  if  Ac  £  6.25 
Ac 


if  1.6  <  Ac  <  6.25 


FrcLr:  ^’FreLr  ^  =  0.329 

FrcLr=3853x  1()8pa 
d  -  Failure  Pressures 

PcLr:=Tl'a0L‘PeL 

|Jrt709^pi 
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4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 
4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


# 


# 


* 


axL:= 


0.207  if  —  £  1242 
t 

(169)  D  M  . 

v  —  if  —  <  1242 


195  + 


••(?)  1 


axG:= 


0.72  if  Arb  >  0.2 

[(3-6  -  5.0-axIj  Arb  +  if  0.06  <  Arl 
if  Afy  <  0.06 


Fxe&:=“*G°6°5E(!}(1*A.b) 


FxeGr=  >•>*  10  p« 


c  -  Inelastic  Buckling  Stresses 


Ac  := 


rxeGr 


n  := 


1  if  Ac  £  0.55 

—  +  0.18^  if  0.55  <  Ac  £1.6 
Ac  ) 

1.31 


1  +  1.15Ac 
1 


if  1.6  <  Ac  <6.25 


if  Ac  >  6.25 


Ac 

*xcGr:=rl  ^xeGr 
FxcGr=3-791*  1()8pa 


< 


<  0.2 
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Buckling  Stresses  With  or  Without  End  Pressure 


if  Mx>  1.56 


\  if  Mx£  1.56 


Ier:=Ir+Arz/ 


V* 


Ar+Le-t 


n  :=6 


Given 
2  ^  n  <  15 

nj  :=  Minimizc^peQ,n) 
nj  =  3.984 


KeG:=(l-0.3k). 


Le‘ 


Af  +  Le*t 


a0G:=°.S 
FreGr:=  a6G' 


PedM-Ko 


Keo 


FreGr=2  04x  10  Pa 


c-  Inelastic  Buckling  Stresses 
^reGr 

AC  :== - 


r\  := 


1  if  Ac  <  0.55 


(—  +  0.181  if  0.55<Ac<n.e 

l  AC  J 


1.31 


1  +  l.lSAc 
1 


if  1.6<  Ac  <  6.25 


if  Ac  ^  6.25 


Ac 


FrcGr:=TlFreGr 

FrcGr=431*  1()8pa 
d  -  Failure  Pressures 

PcGr:=1V“0G-PeG(nl) 


SUMMARY 

4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 

4.1.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


FxeLr=  1-886*  10  Pa 
b  -  Inelastic  Buckling  Stresses 

FxcLr=  3984  x  1()8pa 
4.1.2  -  External  Pressure 

a  -  Elastic  Buckling  Stresses 
FreLr=  1.172x  109Pa 
c  -  Inelastic  Buckling  Stresses 
FrcLr=  3.853  x  108pa 
d  -  Failure  Pressures 


FxeLr=  1-886x  10  Pa 


FxcLr=  3.984x  10  Pa 


FreLr=  1.172x  10  Pa 


FrcLr=3853x  10  Pa 


PcLr  =  4  895x  10  Pa  |^=-70?.97psi 

4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 

4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


FxeGr=l  lx  10  pa 
c  -  Inelastic  Buckling  Stresses 

FxcGr=3-791x  lo8pa 

4.2.2  -  External  Pressure 


FxeGr=l  lx  10  Pa 


FxcGr=3791x  10Pa 


a  -  Buckling  Stresses  With  or  Without  End  Pressure 


FreGr=2  039x  10  Pa 
c  -  Inelastic  Buckling  Stresses 

FrcGr=4-314x  l°8pa 
d  -  Failure  Pressures 


FreGr=2  039x  10  Pa 


=  4.314x  10  Pa 


pcGr=7.39x  10°  Pa 


$cGr- 1071-8Psi 
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Classification  Society  Solution 
(API  Bulletin  2U) 

NAVSEA  Test  Cylinder  2.c 


# 


# 


General  Defintions 

ksi  :=  6.894757106Pa  rtog  :=  64.0—  kip  :=  lOOdbf 

ft3 


E  :=  30000 ksi 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison’s  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 


bload  :=  1 


kpa  := 


kip 

in 

kg 

2 


sec  m 


diameter  to  outside  of  shell  D0  =  38.102in 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


^  =  0.127^1 
hw  =  2.01in 


bf  =  1.552in 
tf  =  0.305in 
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radius  to  outside  of  shell 


R  :=  Rn  —  radius  to  centerline  of  shell 
^  2 


+  tfbr  h.„  +  — 


—  +  if  of  ^nw 

Vhw+tfbf 


2  J  dist  from  shell  to 


centroid 


c,  =  1.757in 


c2 hw  “  C1 


dist  from  centroid  to  end  of  flange 


C2  =  0.253in 


zr  :=(ci  +  -\(-l)  distance  from  centerline  of  shell  to  centroid 

V  2/  of  ring  stiffener  (positive  outward )  Zr  =  — 1.925ii 


\  of  nng  stiffener  (positive  outward ) 

R,.  :=  R  +  Zr  radius  to  centroid  of  ring  stiffener 


Rj  =  16.957in 


Ar  :=  +  bf  tf)  cross-sectional  area  of  ring  stiffener  Ar  =  0.729in2 


L  :=  Lj.  -  b 


unsupported  shell  length 


L=3.129in 


1  3  f  1  3  ^ 

*r :=  ^  VNv  +  Vbw ‘^cl  “  ~j  +  ~bf tf  +  bf tf  c2  +  Jj 


Ir  =  0.312in 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


6  :=  v  3-(  1  —  v/—=- 


Le  :=  1.56^R^tf  cosh  (e) — pos(e)^  effective  shell  length,  [P&S  eqn  92] 

V  sinh(0]  +  sin(G)  ) 


Le  =  3.03in 


.  ,  t  lf 

d:=h  +  -  +  — 
2  2 
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Ap  Le-t 


Aw:  Afl  • 


At:=ap  +  Aw+  Afl 


Aw^ 

~J 


+  AffAp 


yP:=2t+d‘(1_C2) 


yf 


ltf+d-C2 


Ie:=Ard  Cj 


Ie  =  1.913in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


*c 


radius  to  centroidal  axis  of  combined  ring 
stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


R^,  =  18.082in 


# 
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Failure  Modes 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 
4.1.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

Mx:=-p= 

yfBA 

D:=2R 


Mx=  1.291 


D  =  37.765in 


c  := 


2.64  if  MX£1.5 


3.13 


K) 


0.42 


if  1.5<  Mv<  15 


1.0  if  Mx^  15 


axL:= 


0.207  if  —  ^  1242 
t 


(169c) 


195+0 


if  —  <  1242 


c  =  2.64 


«xL=  1777 


Cx:=  0.605 


FxeLr:_  axL^x7^ 


for  D/t  >  300 
t 


D 


rxeLr 


=  3.969x  10  Pa 


b  -  Inelastic  Buckling  Stresses 


rxcLra’- 


233F„ 


166+  1 
(o.5Fy)  if  600 


if  —  <  600 


rxcLr 


:=min((F: 


xeLr  rxcLra, 


» 


FxcLr=  U36x  10  Pa 
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Choose: 

k  =  0.5  k  =  0  for  radial  pressure  and  0.5  for  hydrostatic  pressure 


a  -  Elastic  Buckling  Stresses 


-  =  112.062 
t 


Am:=Mx-L17+  1068k 


Cp  := 


(2Am) 

(t) 


Cp:=Cp 


Am=  0.655 


Cp  =  0.012 


Mx=  1.291 


PeL: 


l°s_.E/±V 

ns. _  U; 


V  +w 


v  w 


if  Mx>  1.5  a  Am<  2.5 


if  2.5<  Am<  0.104^ 


e.es&Cp'1-061-^^ 


if  0.208  <  Cp  <  2.85 


Mi) 


if  Cp  >  2.85 


Mx  is  too 
small, 
notin 
range 


"Value  not  calculated 
(Mx  outside  limits) 


a9L  ;=  o.8  imperfection  factor,  normally  0.8  for  fabricated  cyl 


A:=Ar- 


A  =  0.903in2 


section  1 1 .3 


L  :=  1.56a/R4  +  t, 


w 


Le  =  4.062in 


y:= 


1.0  if  Mx<  1.26 


(1.58-0.4^  if  1.26<MX<3.42 


0  if  Mx>  3.42 


e  :=  • 


1  -  0.3k 


1  + 


V* 
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_  “eL'PeL’fy  „ 

FreLr:“  - - %L 

FreLr=,Pa 


c  -  Inelastic  Buckling  Stresses 


1  if  Ac  <0.55 
( 0.45  't 

- +  0.18  if  0.55  <  Ac  <,  1.6 

V  Ac  ) 

- — -  if  1.6  <  Ac  <  6.25 

1+  1.15Ac 

—  if  Ac  >  6.25 
Ac 

FrcLr:=Tl'FreLr 

FrcLr= ' Pa 
d  -  Failure  Pressures 

PcLr:=1la0LPeL 


**Value  not  calculated 
(Mx  outside  limits) 


4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 
4.2.1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 

A  • 

Arb  •  t 


axL:= 


0.207  if  —  £  1242 
t 

(169)  D 

v  if  —  <  1242 


195+0.5- 


(D)  t 

U  J 


axG:= 


0.72  if  Arb  >  0.2 

[(3.6-5.0axI)-Arb+aJ  if0.06<Ar] 
axL  if  Arb  -  0  06 


FxeGr:=<>xG»^&(j}(1  +  Arb) 


FxeGr=2074*  10  Pa 


c  -  Inelastic  Buckling  Stresses 


Ac  :=■ 


rxcGr 


r\  := 


1  if  Ac  <  0.55 

f—  +  0.18!  if  0.55 <  Ac  <1.6 
l  Ac  J 


1.31 


1+  l.lSAc 
1 


if  1.6  <  Ac  <6.25 


if  Ac  £  6.25 


Ac 

FxcGr:=  ’l  FxeGr 


FxcGr=8  482x  10  Pa 


<0.2 
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Buckling  Stresses  With  or  Without  End  Pressure 


Le== 


(l.lVlM  +  tw)  if  Mx>  1.56 
Lj.  if  Mx^  1.56 


*er '  V  +  Ar'^r  * 


2  Le-t  V‘ 
+ 


Ar  +  V* 


12 


n  :=  6 


Given 
2  £  n  <  15 

nj  :=Minimiz^peQ,n) 


nj  =  2 


KeG:=(l-0.3k)- 


V 

Ar  +  Le-t 


a0G:=08 

FreGr:=a0G 


FreGr=  148x  lo9pa 


c  -  Inelastic  Buckling  Stresses 

FreGr 
Ac  := - 

Fy 

1  if  Ac  <  0.55 

f  —  +  0.18)  if  0.55 <  Ac  <  l.( 
l  Ac  J 

1  31 

-  if  1.6  <  Ac  <  6.25 

1  +  l.lSAc 

—  if  Ac  S  6.25 
Ac 

*VcGr:” ^  *reGr 
FrcGr=7.53x  108Pa 

d  -  Failure  Pressures 

PcGr:=Tl'a80!,eC^nl) 

Bllllill 


RESULTS 


4.1  -  Local  Buckling  of  Unstiffened  or  Ring  Stiffened  Cylinders 

4.1 .1  -  Axial  Compression  or  Bending 
a  -  Elastic  Buckling  Stresses 


FxeLr=  3  969x  10  Pa 
b  -  Inelastic  Buckling  Stresses 

FxcLr=  U36x  1()9pa 
4.1.2  -  External  Pressure 

a  -  Elastic  Buckling  Stresses 


=  3.969x  10  Pa 


‘xeLr 


FxcLr=  1136x  10  Pa 


FreLr=*Pa 

c  -  Inelastic  Buckling  Stresses 

FrcLr=lPa 
d  -  Failure  Pressures 

PcLr='Pa 

4.2  -  General  Instability  of  Ring  Stiffened  Cylinders 

4.2.1  -  Axial  Compression  or  Bending 
a  •  Elastic  Buckling  Stresses 

FxeGr=  2074 x  1()9pa  FxeGr=2074x  1()9pa 

c  -  Inelastic  Buckling  Stresses 

FxcGr=8482x  1()8pa  FxcGr=8482x  1()8pa 

4.2.2  -  External  Pressure 

a  -  Buckling  Stresses  With  or  Without  End  Pressure 

FreGr=  1-479x  ^ Pa  FreGr=  1479 x  l«9pa 

c  -  Inelastic  Buckling  Stresses 

FrcGr=  7534x  1()8  Pa  FrcGr=  7534x  10 *  Pa 

d  -  Failure  Pressures 

PcGr=2-609x  10?Pa  |cGr=3784  63psi 


FreLr=,Pa 


FrcLr=,Pa 


>icLr  =  •  psi  **Value  not  calculated 

(Mx  outside  limits) 
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Classification  Society  Solution 
(PNV-RP-C202) _ 


NAVSEA  Test  Cylinder  1  .d 


Buckling  Strength  of  Shells 
Oct  2002 


General  Defintlons 


ksi  :=  6.89475710  Pa 


E  :=  30000 ksi 


v  :=  .3 


_  „lbf 
rtog  :=  64.0 — 

ft3 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison's  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 

D0  :=  2-(R|  +  t)  diameter  to  outside  of  shell 


kip  :=  lOOflbf  bload:=l-; 

kpa  := 


kip 


in 

kg 


2 

sec  -m 


D0=  16.176in 


ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


b:=V 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw  =  0.138in 
hw  =  0.57in 


bf  =  Oin 
tf  =  Oin 


■v3 

R:=Ro-i 


radius  to  outside  of  shell 


radius  to  centerline  of  shell 


R  =  8.047in 


Kv'^w  r 


- +  tfbf  h...  h - , 

2  t  t  ^  w  2  J  dist  from  shell  to 

°i:-  Vhw  +  tfbf  centroid 


C|  =  0.285in 


c2  :-hw_cl 


dist  from  centroid  to  end  of  flange 


C2  =  0.285in 


Zf  :=  cj  +  -  distance  from  centerline  of  shell  to  centroid 
2  of  ring  stiffener  (positive  outward ) 

Rj.  :=  R  +  Zf  radius  to  centroid  of  ring  stiffener 
Aj.  :=  +  bf  tf)  cross-sectional  area  of  ring  stiffener 


Zf  =  0.325in 


Ar  =  0.079in2 


L  :=  Lj.  -  b 


unsupported  shell  length 


L  =  4.128in 


1  3  f  KY  1  3  (  tfV 

V-  22 +  ^C1  Y)  +  ~\2 +  c2  + 


^  A 

Ir  =  2.13x  10  in  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


0  :=  v3  \l  -  v2) — — 


Le  :=  1  -56sJlTt  ( effective  shell  length,  [P&S  eqn  92] 
V  sinh(e)  +  sin(e)  J 


Le  =  1.255in 


t  *f 
d:=hw+-  +  - 
2  2 
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Afl  :=  tf  bf 


At  :=  Ap  +  Aw  +  Afl 


Ie:=Axa  Cj 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


Ie  =  0.007in4 


radius  to  centroidal  axis  of  combined  ring 
stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


Rp  =  8.181in 
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Failure  Modes 

3.5  Ring  stiffened  shells 


3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling 


V  :=2  p  :=.6 
$5=1.04^ 


Buckling  coefficients  for 
hydrostatic  pressure 
(Table  3.4-1) 


%h:=fE 


vsV 

<  V  ) 

-viz'2 

Elastic  buckling  strength  of 

Km 

unstiffened  shell 

(3.4.1) 

fEa:=fE 

(3.4.2) 


%  =  2.555x  10  Pa 


a 


P  := 


\  • 

1.56VR4 


sinh(p)-cos(p)  +  cosh(p)sin(p) 
sinh(2-p)  +  sin(2-p) 


m 


0  if^<0 


(2.2.10) 
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Psd  SB  := 


press  <—  lOOpsi 
limit  4-  5psi 
test  <-  O  psi 
convert  4-  1-psi 
j<-0 

while  j  <  40 


-press  R 

ua_sd 

2t 

-press  *R 

°h_sd 

t 

( 

2 

aj_sd  v 

^°a_sd 

F„, 

A's_square 

y 

<— 

aj_sd 

F., 

n  . _ 

y 

!  i  v>l  e 

°i 

a  +  1 


CTa_sd‘CTh_sd  +  CTh_sd  ) 

^  “^a  sd  “^h  sd  ^ 
- = —  + 


fEh  J 


f' 


1  +  ^s_square 
s  '  i/^ssquare 


X, 


1.15  if  Xg  <  0.5 


(0.85+  0.&Xsj  if  0.5  <  Xg  ^  1.0 


1.45  if  Xg  >  1.0 


fksd  4 - 


ks 

ym 

break  if  aj_sd  >  f^ 

press  <-  press  +  2*psi 
j<-j+i 

press 

out  4 - 

u  convert 

OUt  j  4—  j 


OUt2  4-  Xg 


out 


(2.2.2) 


(2.2.14) 


(3.2.3) 

(3.2.2) 

(3.2.1) 


(3.1.3) 


(3.1.2) 

(3.1.1) 
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Psd  SB 


(  166  ^ 

33 

1.934,/ 


Psd_L  :=  Psd_SB0-  !psi  “  2Psi 


Maximum  pressure  to  still  meet  stability 
requirement  of  eqn.  3.1.1. 

Prevent  shell  shell  buckling  (Local  Buckling) 


3.5.2  Panel  ring  buckling 


3.5.2. 1  Cross  sectional  area  check 


The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


Rt 


vReq := 


—  +  0.06^  Lj.  t  (3.5.1) 


VZL 


required  area 


actual  area 
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3.5.2.7  Refined  calculation  of  l_h  for  external  pressure 

Method  for  calculating  the  the  capacity  of  the  ring  frame 


If  a  ring  stiffened  cylinder,  or  a  part  of  a  ring  stiffened  cylinder,  is 
effectively  supported  at  the  ends,  the  following  procedure  may  be  used 
to  calculate  required  moment  of  inertia. 

Moment  of  inertia  for  the  combined  plate/stiffener 
previously  calculated  [Hughes  eqn  8.3.6] 


V=ie 


Ih  =  0.00658in4 


fj- :~Fy 
zt:=yf 

:=  0.005R 

r  l.56yJvTt') 


characteristic  material  strength  (yield  strength) 


*eo 


min  * 


^  ) 


*eo  :=  mic(leo  min)  equivalent  length 


leo  =  1259m 


a 


(3.5.24) 


D0 
r  :=  — 
r  2 

rf:=R+  2  +  hw+tf 
\  :=Psd_SB2 


radius  of  shell  measured  to  ring  flange 
(external  stiffeners  +  +  +  /  internal  -  -  -  ) 


ym- 


1.15  if  xs  <  0.5 


material  factor  from  sec  3.1 
(0.85+0.60A.S)  if  0.5<XS  ^  1.0  (3 A  ^ 


1.45  if  Xs  >  1.0 


Ym=145 
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Square  :"Ar+leo.t 

12  (l  -  v2)-Ih 


aB;=- 


ci:= 


2(1  +  <xb)/ 


'h_squarc  ~  0.036in 


aB  =  31.671 


1  +  a 


0.27-Zt  aB  ^ 

1  + . . C,  =  185.478 


J1  +  aB  1  +  aB )  1 


C2:=2^1  +  0.27Zl 


Cj  =  28.342 


H:=- 


Zt^ o*Th 
*h_square  **Vo 


j_  92^  1 


Cl)  l-l 
2 


f  :=Cl-f  E  T/-i 
"  112.(l-v2)U 


hj 


H  =  1.917 


%=  1.81  lx  107psi 


Xl:=  I7 
% 


2lj  =  0.066 


1+  M  +  J 


(l  +  p  +  X,2)  -  4-Xj2 


2Xj 


1  + 


Psd_GI:=0/75 - trf 


M 


fjc=  2.74  x  10*  psi 


(3.5.27) 

(3.5.23) 

(3.5.21) 

(3.5.28) 

(3.5.25) 

(3.5.20) 

(3.5.16) 

(3.5.15) 

(3.5.14) 


maximum  allowed  external  pressure  to  prevent  par 
ring  buckling  (General  Instability) 


SUMMARY 


3.5  Ring  stiffened  shells 

3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling  (Elastic  local  buckling  of  unstiffened  shell) 


3.5.2  Panel  ring  buckling 

3.5.2.1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area  actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure  (General  Instability) 
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Classification  Society  Solution 
(DNV-RP-C202) 

Buckling  Strength  of  Shells 
Oct  2002 

NAVSEA  Test  Cylinder  1  .f 


General  Defintions 

lbf 

ksi  :=  6.894757 106Pa  rtog  :=  64.0 —  kip:=100(lbf  bload  := 

ft3 


E  :=  30000 ksi 


v  :=  .3 


Young's  Modulus  of  Elasticity  Kpa  • 

Yield  Strength 

Poison’s  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 

diameter  to  outside  of  shell  D0  =  34.92in 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


in 


kg 


2 

sec  -m 


Ring  Stiffener  Dimensions 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 


flange  thickness  of  ring  stiffener 


tw  =  0.198in 
hw=0.762in 


bf  =  0.763in 
tf  =  0.263in 


R»R0-i 


radius  to  outside  of  shell 


radius  to  centerline  of  shell 


R  =  17.328in 


,  ft 

- +  tfbf  h...  +  —  , 

2  i  t  ^  w  2  J  dist  from  shell  to 

Vhw+tfbf  centroid 


Cj  =  0.674in 


c2:-hw-ci 


dist  from  centroid  to  end  of  flange 


Zr  :=  fcj  +  -Y(-i)distance  from  centerline  of  shell  to  centroid 
V  2  J  of  ring  stiffener  (positive  outward ) 

Rj  :=  R  +  Zt  radius  to  centroid  of  ring  stiffener 
Aj.  :=  (twhw  +  bftf)  cross-sectional  area  of  ring  stiffener 


C2  =  0.088in 


Zr  =  -0.805in 


Af  =  0.352in 


L^L^-b 


unsupported  shell  length 


L  =  2.468in 


1  3  (  1  3  f  *fV 

V:=  ^Vhw  +  W^cl  “  —J  +  ^*bf  *f  +  bf ‘f  |^c2  +  “j 


I,.  =  0.03  tin 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener  • 


0  :=v3-\l  -  v 


2  L 


:=  1.56^R^t  f -sh-|9^ — ^  effective  shell  length,  [P&S  eqn  92] 
V  sinh(0)  +  sin(e)  ) 


Le  =  2.41  in 


t  *f 
d:=hw+-  +  - 
w  2  2 
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Ap  Le-t 


Aw 


Afl :=  tf  bf 


AT:=Ap  + Aw+ Afl 


Cl:: 


Aw^ 

-Tj  +  AffAp 


AW  ^ 

- +  An  * 

l  2  V 


Aj 


yP 


1 

— 1+ 

2 


d-fl-Cz) 


yf 


ltf+d-02 


Ie  :=  Aj-d2^! 
Ie  =  0.182in4 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


*c 


radius  to  centroidal  axis  of  combined  ring  -  i7.041in 

stiffener  and  effective  width  of  shell  c 

(outside  stiffener  -/+.  inside  +1- 1 
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Failure  Modes 


3.5  Ring  stiffened  shells 
3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling 


Buckling  coefficients  for 
hydrostatic  pressure 
(Table  3.4-1) 


(3.4.2) 


9sb'*2,e 

12-(l  -  v2)  1m 


Elastic  buckling  strength  of 
unstiffened  shell 
(3.4.1) 


fEh:=fE  %a;=fE 


%  =  3.893  x  109Pa 


>eo  :=  1.56VR1 


a 


P:= 


\  • 

l.5$y[RA 


sinh(p)cos(p)  +  cosh(p)sin(p) 
sinh(2-p)  +  sin(2-p) 


m 


0  if  4m<0 


(2.2.10) 
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Psd  SB  := 


press  <-  1600psi 
limits  5  psi 
test  <-  O  psi 
convert  l  psi 
j<-0 

while  j  <  100 


-^press  *R 

CTa_sd 

2-t 

-press  *R 

CTh_sd 

t 

I 

f  2 

aj_sd  \ 

^CTa_sd  “ 

Fv 

<_  y  . 

ssquare 

CTj_sd 

F 

n 

y 

1- 


“('-tK- 


a  +  1 


.  2 


aa_sd*CTh_sd  +  ah_sd 
~CTa_sd  ~~<Th_sd  ^ 

fEa  %h  J 


J 


1  +  X, 


ssquare 


^s  -y/^s  square 


1.15  if  Xs  <  0.5 


^0  85  +  0.6-A.sj  if  0.5  <  A,s  <  1.0 


1.45  if  Xs  >  1.0 


*k s 


*ksd 

Ym 

break  if  Gj_sd  >  f^ 

press  <-  press  +  2  psi 
j<- j+  1 

press 


0UV 


out. 


out. 


convert 


J 


out 


(2.2.2) 

(2.2.14) 

(3.2.3) 

(3.2.2) 

(3.2.1) 

(3.1.3) 

(3.1.2) 

(3.1.1) 
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Psd_L:=Psd_SB0' 


lpsi  -  2psi 

Maximum  pressure  to  still  meet  stability 
requirement  of  eqn.  3.1.1. 

Prevent  shell  shell  buckling  (Local  Buckling) 


3.5.2  Panel  ring  buckling 


3.5.2. 1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area  actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure 

Method  for  calculating  the  the  capacity  of  the  ring  frame 


If  a  ring  stiffened  cylinder,  or  a  part  of  a  ring  stiffened  cylinder,  is 
effectively  supported  at  the  ends,  the  following  procedure  may  be  used 
to  calculate  required  moment  of  inertia. 


Moment  of  inertia  for  the  combined  plate/stiffener 
previously  calculated  [Hughes  eqn  8.3.6] 


,  4 


Ih  =  0.18238in 


V=Ie 

fj.  :=  Fy  characteristic  material  strength  (yield  strength) 

zt:=yf 
<^:=  0.005 R 

/  L56^s7t\ 


eo  mm  ■ 


V 


**  ) 

min)  equivalent  length 


leo  =  2.666in 


a  := 


leo'* 


Dn 


r,:=- 


rf:=R“  2  -hW"-tf 


h  :=Psd_SB2 


(3.5.24) 


radius  of  shell  measured  to  ring  flange 
(external  stiffeners  +  +  +  /  internal  —  ) 


7m:= 


1.15  if  <  0.5 


(0.85  +  0.60 if  0.5  <\<  1.0 
1.45 

Ym=  U93 


material  factor  from  sec  3.1 
(3.1.3) 


1.45  if  Xs>1.0 
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^square  :"Ar+Ieo.t 
12  (l  -  v2)  Ih 


aB:=- 


Ci:= 


2-(l  +  aB) 
1  +  a 


2 

'hsquarc  =  0-1 73  in 


aB  =  41.065 


f  I  0.27-Zl  aB  ^ 

1  +  — . . . —  C,  =  172.677 

U  >/1  +  aB  1  +  aB  J 


Cj  :=  2-^1  +  0.27Zl 


20.13 


H  :=  7 


^square  R1eo 


j  1 


l  cUi_* 


p  =  0.422 


%:=Cl— fa 

12  (1  -  v2)  l 


fB  =  2.923  x  107  psi 


^JI 


=  0.058 


*  +  P+>-l  -J 


^1  +  p  +  A.j  ^  —  4-X.j 


2-X, 


fk  =  6.924x  104psi 


1  + 


Psd_GI:-°-75 - trf  ~f  n 

7m  [  ,  v  ] 


*2H 


(3.5.27) 


(3.5.23) 


(3.5.21) 


(3.5.28) 


(3.5.25) 


(3.5.20) 


(3.5.16) 


(3.5.15) 


(3.5.14) 


maximum  allowed  external  pressure  to  prevent  par 
ring  buckling  (General  Instability) 


SUMMARY 

3.5  Ring  stiffened  shells 

3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling  (Elastic  local  buckling  of  unstiffened  shell) 


3.5.2  Panel  ring  buckling 


3.5.2. 1  Cross  sectional  area  check 


The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area 


actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure  (General  Instability) 
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Page  Intentionally  Left  Blank 


{Classification  Society  Solution 
(DNV-RP-C202) _ 


NAVSEA  Test  Cylinder  2.a 


Buckling  Strength  of  Shells 
Oct  2002 


General  Defintions 


ksi:=  6.89475710  Pa 


rtog  :=  64.0- 


lbf 


v  :=  .3 


ft 


E  :=  30000 ksi  Young's  Modulus  of  Elasticity 


kip  :=  lOOdbf  bload  :=  1  - 
kpa  := 


kip 


in 

kg 


2 

sec  -m 


Poison’s  ratio  for  Fe/Steel 
thickness  of  shell 


Inner  radius  of  cylinder 
Dq  :=  2-^  +  t)  diameter  to  outside  of  shell 


Dq  =  16.922in 


ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


Ring  Stiffener  Dimensions 


b  :=t 


w 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 

flange  thickness  of  ring  stiffener 


tw  =  0.044in 


hw  =  0.454in 


bf  =  0.399in 


tf  =  0.078in 
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*0»- 


RWR(,-- 


*w'b’ 


w 


radius  to  outside  of  shell 
radius  to  centerline  of  shell 


c,  :=• 


f 

+  tfbf|hw+  — 


A 

) 


Vhw  +  tfbf 


dist  from  shell  to 
centroid 


R  =  8.418in 


Cj  =  0.389in 


c2;=hw-ci 


dist  from  centroid  to  end  of  flange 


C2  =  0.065in 


Rj  * —  R  +  Zj. 


distance  from  centerline  of  shell  to  centroid 
of  ring  stiffener  (positive  outward ) 

radius  to  centroid  of  ring  stiffener 


Aj.  :=  (twhw  +  bftf)  cross-sectional  area  of  ring  stiffener 


Zr  =  0.432in 


Ar  =  0.051m 


L  :=  Lj.  -  b 


unsupported  shell  length 


L=  1.322in 


.1,3 
j2**w*”w  + 


f 


bw^  1  3 

Ci - 1  H - -bftf  +  bftf 

V  1  2  )  12  1  f  f  f 


Co  H - 


jfV 

2) 


^  A 

\  =  1.219x  10  in  moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


0:=j3-(l-v2)-pr 

VrI 

Le:=1.56^/^#l£2#"i 

V  sinh(o)  +  sin(o)  ) 


effective  shell  length,  [P&S  eqn  92] 


Le  =  1.221in 
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Ap  Le-t  Aw A^ .  tp  bp  A-p .  Ap  +  Aw  +  Ajj 


AT  Aw^I 


w  I 


Ci:= 


3  4 


-j  +  ftflAp 


C2:=- 


(T^) 


yp:=lt  +  d.(l-C2) 


yf:=  2tf+dC2 


re :=  AT'd  C1 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


Ie  =  0.008in 


RC:=R--  +  yp 


radius  to  centroida!  axis  of  combined  ring 
stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+.  inside  +/- ) 


Rg  =  8.559in 
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Failure  Modes 

3.5  Ring  stiffened  shells 
3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling 


Buckling  coefficients  for 
hydrostatic  pressure 
(Table  3.4-1) 


(3.4.2) 


,  Csfi-n  E  (  t  V 

— / - c-  — 


Elastic  buckling  strength  of 


i  „2)  I  4  J  unstiffened  shell 
(3.4.1) 


%=  1.643x  10  Pa 


%h:=fE 


fEa:=fE 


^■.=  1.56^1 
A, 


^eo’* 


P:= 


1.56^/iTt 


_  __2  sinh(p)  cos(p)  +  cosh(p)sin(p) 
sinh(2-p)  +  sin(2-p) 


0  if  ^<0 


(2.2.10) 
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Psd  SB  := 


press  <—  600 psi 
limit  5*  psi 
test  <-  O  psi 
convert  <-  1-psi 

j<-0 

while  j  <  40 

-press  R 
2-t 


Ta  sd 


7h  sd 


Tj_sd 


-press  -R 


1- 


a  +  1 


(a< 


a  sd  CTa  sd*ah  sd  +  ah 


A 

i  sd  J 


ssquare 


aj__sd 


”aa_sd  "^h  sd  ^ 
V  %a 


%h  J 


lks 


J 


i  + 1 


s_square 


*■»  <- 


1.15  if  Xg  <  0.5 


(0.85+0.6-A.S)  if  0.5  :£  £  1.0 


1.45  if  Xs  >  1.0 


rks 


rksd 


break  if  aj  sd  > 
press  <-  press  +  2*psi 

j<-  j+ 1 

press 

out  < - 

u  .  convert 

outj  <- j 

out2  h 


out 


(2.2.2) 

(2.2.14) 

(3.2.3) 

(3.2.2) 

(3.2.1) 

(3.1.3) 

(3.1.2) 

(3.1.1) 
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Psd  SB  = 


39 

K0.1\l) 


Psd_L:=Psd_SB01Psi-2Psi 


Maximum  pressure  to  still  meet  stability 
requirement  of  eqn.  3.1.1. 

Prevent  shell  shell  buckling  (Local  Buckling) 


3.5.2  Panel  ring  buckling 

3.5.2. 1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area 


actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure 


Method  for  calculating  the  the  capacity  of  the  ring  frame 


If  a  ring  stiffened  cylinder,  or  a  part  of  a  ring  stiffened  cylinder,  is 
effectively  supported  at  the  ends,  the  following  procedure  may  be  used 
to  calculate  required  moment  of  inertia. 


Moment  of  inertia  for  the  combined  plate/stiffener 
previously  calculated  [Hughes  eqn  8.3.6] 


V=Ie 

<r:=Fy 

zt:=yf 

:=  0.005R 


Ih  =  0.00773in4 

characteristic  material  strength  (yield  strength) 


Jeo 


min  * 


^  ) 
Igp  :=  ™in(leo_mm) 


equivalent  length 


leo  =  l-326in 


a  := 


Wo'1 


(3.5.24) 


fr'  2 

t  radius  of  shell  measured  to  ring  flange 

rf  :=  R  +  —  +  hw  +  tf 

2  (external  stiffeners  +  +  +  /  internal  —  ) 

V=Psd_SB2 


Ym:= 


1.15  if  Xs  <  0.5 


(0.85  +  0.60 if  0.5  <X&<  1.0 


material  factor  from  sec  3.1 
(3.1.3) 


Hi  square 


^square 


=  0.047in2 


aB  =  97.872 


2-(l  +  aB) 

0.27-Zt 

1  +  L 

1  +  a 

y  ^ 

aB  ^ 

-  Ci  =  126.503 

l  +  aBj 


1  +  0.27- ZT 


C2  =  10.506 


M  := 


’hsquare  '^'eo 


,  1 


CU  i-i 


(i  =  0.462 


f  —  r  »2~e  (  <  \ 


fB=  3.408x  107  psi 


Xl:= 


Xj  =  0.044 


%  =  (*r) 


1  +  p  +  M  -J 


(l  +  p  +  A.,2)  -  4-Xj2 


2-Xj 


fjc  =  4.479x  10  psi 


1  + 


Psd  GI:_0-75 - t  rf 


ym 


H) 


(3.5.27) 


(3.5.23) 


(3.5.21) 


(3.5.28) 


(3.5.25) 


(3.5.20) 


(3.5.16) 


(3.5.15) 


(3.5.14) 


maximum  allowed  external  pressure  to  prevent  par 
ring  buckling  (General  Instability) 


SUMMARY 

3.5  Ring  stiffened  shells 

3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling  (Elastic  local  buckling  of  unstiffened  shell) 


3.5.2  Panel  ring  buckling 


3.5.2. 1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area  actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure  (General  Instability) 
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Buckling  Strength  of  Shells 
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NAVSEA  Test  Cylinder  2.c 


General  Defintions 


ksi:=  6.894757 10  Pa 


rtog  :=  64.0 


lbf 


E:=  30000 ksi 


v  :=  .3 


D0»2(Ri+l) 


r 


kip  :=  100 Clbf  bload  :=  1  — 

kpa :: 


kip 


Young's  Modulus  of  Elasticity 
Yield  Strength 

Poison's  ratio  for  Fe/Steel 
thickness  of  shell 

Inner  radius  of  cylinder 

diameter  to  outside  of  shell 

ring  spacing  (frame  center  to  frame  center) 

length  of  cylinder  between  bulkheads  or  lines  of  support 


in 

kg 


2 

sec  m 


Dq  — -  38.102in 


Ring  Stiffener  Dimensions 


b :  tyj 


thickness  of  web  of  ring  stiffener 

height  of  web  of  ring  stiffener 

faying  width  of  stiffener  (from  P&S 
for  I  beam  stiffener) 

breadth  of  flange  of  ring  stiffener 


tw  =  0.127in 
hw  =  2.01in 


bf=  1.552in 


flange  thickness  of  ring  stiffener 


tf  =  0.305in 


D  °o 
*°:=T 

R:=Ro-i 


radius  to  outside  of  shell 


radius  to  centeriine  of  shell 


R  =  18.883in 


+  tft)f  hw+  — 


2  )  dist  from  shell  to 


centroid 


Cj  =  1.757in 


c2:-hw-cl 


dist  from  centroid  to  end  of  flange 


zr :=  f  cl +  -l  (-l)distance  from  centerline  of  shell  to  centroid 
V  2)  of  ring  stiffener  (bositive  outward^ 

Rj  :=  R  +  Zr  radius  to  centroid  of  ring  stiffener 
Aj.  :=  (twhw  +  bftf)  cross-sectional  area  of  ring  stiffener 


C2  =  0.253in 


Zr  =  -J.925in 


Ar=0.729in2 


L^Lj-b 


unsupported  shell  length 


L=3.129in 


1  3  (  1  3  (  tf7 

*r :=  '^'Vhw  +  Vhw'^cl  ~~ J  +  72  tf  +  +  ~2) 


If—  0.312in 


moment  of  inertia  of  ring  stiffener  about  its  centroidal  axis 


Moment  of  Inertia  Calculations  for  combined  plate  and  stiffener 


e:=-J3.\l-v 


2  L 


L  :=  1 .56^1 U-f  cosh  (9)  cos(e)7  effective  shell  length,  [P&S  eqn  92] 
V  sinh(9)  +  sin(0j  J 


Le  =  3.03  in 


t  ‘f 
d:=hw+-  +  - 
2  2 
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Ap  Le-t 


Aw  :=tw‘hw  Afl:=tfbf 


+  AflAp 


y  »l'  +  d.(l-C,) 


AT:=Ap  + Aw+ Afl 


yf:=  2'tf+  dC2 


Ie:=ATd  Cj 


Ie  =  1.913in 


moment  of  inertia  for  combined  plate/stiffener 
[Hughes  eqn  8.3.6] 


radius  to  centroidal  axis  of  combined  ring 
stiffener  and  effective  width  of  shell 
(outside  stiffener  -/+■  inside  +/- ) 


=  18.082in 
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Failure  Modes 

3.5  Ring  stiffened  shells 
3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling 

V  :=  2  p  :=  .6 

4:=1.04^zi 
cSB:='t'  | 


Rt 


fEh:=fE 


fp-sf 

U ) 

r(rV 

Elastic  buckling  strength  of 

unstiffened  shell 

(3.4.1) 

%a:=fE 

Buckling  coefficients  for 
hydrostatic  pressure 
(Table  3.4-1) 


(3.4.2) 


leo^l.56  V^t 

Ar 


a  :=  * 


W=2‘ 

m := 


P:= 


1.56  ViTt 


sinh(p)-cos(p)  +  cosh(p)sin(p) 
sinh(2p)  +  sin(2p) 

0  if  ^<0 


fg  =  4.304 x  10  Pa 


(2.2.10) 
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Psd  SB 


press  <—  3100psi 
limit  <-  5  psi 
test  <r-  Opsi 
convert  <-  l*psi 

j<-0 

while  j  <  100 

-press  R 
2*t 


sd 


ah  sd 


<- 


-press  R 
t 


1- 


“  I  1_  2/^-m 


a  +  1 


(  2- 
Vaa_sd  aa_s 


2 ) 

i_sd  J 


aj_sd  ^CTa_sd  CTa_sd'°rh_sd  +  ah 

~  -  ' 


^s_square 


G: 


j_sd 


aa_sd  ""^hsd 
v  %a 


%h  ) 


fks 


J 


1  +  Xt 


s_square 


fit s  square 


1.15  if  Xs  <  0.5 


(0.85+0.6Xs)  if  0.5  <Xg  £  1.0 


1.45  if  \  >  1.0 


rks 


rksd 


Ym 

break  if  CTj  sd  >  f^d 

press  f-  press  +  2  psi 
j<- j+ 1 

press 

out  < - 

u  convert 

outj  <- j 

OUt2  4-  Xs 


out 


(2.2.2) 

(2.2.14) 

(3.2.3) 

(3.2.2) 

(3.2.1) 

(3.1.3) 

(3.1.2) 

(3.1.1) 
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3.122x  10 


Psd_SB=  11  I 
^  0.693  ) 


Psd_L  :=  Psd_SB01Psi  “  2Psi 

_  Maximum  pressure  to  still  meet  stability 

requirement  of  eqn.  3.1.1. 

Prevent  shell  shell  buckling  (Local  Buckling) 


3.5.2  Panel  ring  buckling 

3.5.2. 1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


ZL  := - V 1  -v 

L  R-t 


+  0.06^-Iyt 
) 


required  area 


(3.5.1) 


1  ..  ■  * 

WtfUMWa 

■  -it  '  ’■ 

actual  area 
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3.5.2.7  Refined  calculation  of  l_h  for  external  pressure 

Method  for  calculating  the  the  capacity  of  the  ring  frame 


If  a  ring  stiffened  cylinder,  or  a  part  of  a  ring  stiffened  cylinder,  is 
effectively  supported  at  the  ends,  the  following  procedure  may  be  used 
to  calculate  required  moment  of  inertia. 


Moment  of  inertia  for  the  combined  plate/stiffener 
previously  calculated  [Hughes  eqn  8.3.6] 


V=Ie 

*r:=Fy 

zt:=yf 
:=  0.005R 

Ji 

*eo_min 

V 


Ih  =  1.91274in4 

characteristic  material  strength  (yield  strength) 


.56 

**  j 


igo :  n^leo  min) 


equivalent  length  ^  =  3.256in 


^eo’t 

...  D0 
2 

rf  :=  R  -  —  -  hw  -  tf 


(3.5.24) 


radius  of  shell  measured  to  ring  flange 
(external  stiffeners  +  +  +  /  internal  —  ) 


V=Psd_SB2 


1.15  if  \  <  0.5 

(o.85  +  0.6(US)  if  0.5  <XS  <  1.0 


material  factor  from  sec  3.1 
(3.1.3) 


1.45  if  K  >  1.0 


Ar+,eo‘ 


*h_square 


=  1.048in2 


aB=  167.612 


Ci:= 


2-(l  +  aB) 
l  +  o 


0.27-Zl 
>/1  +  aB 


aB  ^  • 

-  C,  =  1.121x10' 

l  +  aBj 


Cj  :=  2-jr+03JZ^ 


C2  =  46.53 


Wt  (x  <g\  1 

*h_square  'R1eo  ^  Cl)  1  _1 

2 


H  =  0.149 


%:=cr 


g2-E  t  ^ 
12  (1  -v2) 


fg  =  8.869x  107  psi 


=  0.042 


(3.5.27) 


(3.5.23) 


(3.5.21) 


(3.5.28) 


(3.5.25) 


(3.5.20) 


(3.5.16) 


(3.5.15) 


(3.5.14) 


maximum  allowed  external  pressure  to  prevent  par 
ring  buckling  (General  Instability) 


SUMMARY 


3.5  Ring  stiffened  shells 

3.5.1  General 

The  buckling  modes  to  be  checked  are: 

a)  Shell  buckling,  see  Section  3.4.2. 

b)  Panel  ring  buckling,  see  Section  3.5.2. 

c)  Column  buckling,  see  Section  3.8.  -  (not  applicable  for  these  cylinders) 


3.4.2  Shell  buckling  (Elastic  local  buckling  of  unstiffened  shell) 


3.5.2  Panel  ring  buckling 


3.5.2. 1  Cross  sectional  area  check 

The  cross  sectional  area  of  a  ring  frame  (exclusive  of  effective  shell 
plate  flange)  should  not  be  less  than  Areq,  which  is  defined  by: 


required  area  actual  area 


3.5.2.7  Refined  calculation  of  l_h  for  external  pressure  (General  Instability) 
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